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"^ \ Abstract 

i-^ . The aim of these notes is to present in a comprehensive and relatively self-contained way 

some recent developments in the mathematical analysis of two-dimensional viscous flows. 
We consider the incompressible Navier-Stokes equations in the whole plane M. 2 , and assume 
that the initial vorticity is a finite measure. This general setting includes vortex patches, 
vortex sheets, and point vortices. We first prove the existence of a unique global solution, 
for any value of the viscosity parameter, and we investigate its long-time behavior. We 

\J2 . next consider the particular case where the initial flow is a finite collection of point vortices. 

In that situation, we show that the solution behaves, in the vanishing viscosity limit, as a 
superposition of Oseen vortices whose centers evolve according to the Helmholtz-Kirchhoff 
point vortex system. The proof requires a careful stability analysis of the Oseen vortices in 
the large Reynolds number regime, as well as a precise computation of the deformations of 

^— ^ ' the vortex cores due to mutual interactions. 
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1 Introduction 

Although real flows are always three-dimensional, it sometimes happens that the motion of 
a fluid is essentially planar in the sense that the fluid velocity in some distinguished spatial 
direction is negligible compared to the velocity in the orthogonal plane. This situation often 
occurs for fluids in thin layers, or for rapidly rotating fluids where the Coriolis force strongly 
penalizes displacements along the axis of rotation. Typical examples are geophysical flows, for 
which the geometry of the domain (the atmosphere or the ocean) and the effect of the Earth's 
rotation concur to make a two-dimensional approximation accurate and efficient [12] . 

From a mathematical point of view, planar flows are substantially easier to study than 
three-dimensional ones. For instance, it is known since the pioneering work of J. Leray [46(. I47j 
that the two-dimensional incompressible Navier-Stokes equations are globally well-posed in the 
energy space, whereas global well-posedness is still an open problem in the three-dimensional 
case, no matter which function space is used |75j . The situation is essentially the same for the 
incompressible Euler equations, which describe the motion of inviscid fluids [52} 113]. However, 
having global solutions at hand does not mean that we fully understand the dynamics of the 
system. As a matter of fact, we are not able to establish on a rigorous basis the phenomenological 
laws of two-dimensional freely decaying turbulence, and the stability properties of 2D boundary 
layers in the high Reynolds number regime are not fully understood. 

In these notes, we consider the idealized situation of an incompressible viscous fluid filling the 
whole plane M 2 and evolving freely without exterior forcing. Following the approach of Helmholtz 
[36] . we use the vorticity formulation of the problem, which is more appropriate to investigate 
the qualitative behavior of the solutions. Our goal is to understand the stability properties 
and the interactions of localized vortical structures at high Reynolds numbers. This question 
is important because carefully controlled experiments [171 168] . as well as numerical simulations 
of two-dimensional freely decaying turbulence [491 [5], suggest that vortex interactions, and 
especially vortex mergers, play a crucial role in the long-time dynamics of viscous planar flows, 
and are responsible in particular for the inverse energy cascade. Although nonperturbative 
phenomena such as vortex mergers may be very hard to describe mathematically [59] 145] , we 
shall see that vortex interactions can be rigorously studied in the asymptotic regime where the 
distances between the vortex centers are much larger than the typical size of the vortex cores. 

We begin our analysis of the two-dimensional viscous vorticity equation in Section [2j We 
first recall standard estimates for the two-dimensional Biot-Savart law, which expresses the 



velocity field in terms of the vorticity distribution, and we enumerate in Section [2. II some general 
properties of the vorticity equation, including conservation laws, Lyapunov functions, and scaling 
invariance. In Section 12.21 we introduce the space of finite measures, which allows to consider 
nonsmooth flows such as vortex patches, vortex sheets, or point vortices. It is a remarkable fact 
that the vorticity equation is globally well-posed in such a large function space, for any value of 
the viscosity parameter, see [344125] and Theorem 12.81 below. Although a complete proof of that 
result is beyond the scope of these notes, we show in Section [231 that the classical approach of 
Fujita and Kato [24J applies to our problem and yields global existence and uniqueness of the 
solution provided that the atomic part of the initial vorticity distribution is sufficiently small. 
For larger initial data, existence of a global solution can be proved by an approximation scheme 
[16} [Ml S3] but additional arguments are needed to establish uniqueness [25] . 

Section [3] collects a few results which describe the behavior of global solutions of the vorticity 
equation in L 1 (M 2 ). In particular we prove convergence as t — > oo to a family of self-similar 
solutions called Lamb-Oseen vortices, see [31] and Theorem l3.1l below. To establish these results, 
we use accurate estimates on the fundamental solution of convection-diffusion equations, which 
were obtained by Osada [63] and are reproduced in Section 13.11 Another fundamental tool is a 
transformation into self-similar variables, which compactifies the trajectories of the system and 
allows to consider w-limit sets, see Section [3.21 Using a pair of Lyapunov functions, one of which 
is only defined for positive solutions, we establish a "Liouville Theorem", which characterizes 
all complete trajectories of the vorticity equation in L 1 (R 2 ) that are relatively compact in the 
self-similar variables. The conclusion is that all these trajectories necessarily coincide with 
Oseen vortices, see Proposition 13.41 In Section 13.31 we show that Liouville's theorem implies 
Theorem 13. 1[ and proves at the same time that the vorticity equation has a unique solution when 
the initial flow is a point vortex of arbitrary circulation. This is an important particular case of 
Theorem 12. 8[ which cannot be established by a standard application of Gronwall's lemma. 

In Section 0] we investigate in some detail the stability properties of the Oseen vortices, 
which are steady states of the vorticity equation in the self-similar variables. We introduce in 
Section 14.11 an appropriate weighted space for the admissible perturbations, and we show that 
the linearized operator at Oseen's vortex has a remarkable structure in that space : it is the sum 
of a self-adjoint operator, which is essentially the harmonic oscillator, and a skew-symmetric 
relatively compact perturbation, which is multiplied by the circulation of the vortex. This 
structure almost immediately implies that Oseen vortices are stable equilibria of the rescaled 
vorticity equation, and that the size of the local basin of attraction is uniform in the circulation 
parameter. This is in sharp contrast with many classical examples in fluid mechanics, such as 
the Poiseuille flow or the Couette- Taylor flow, for which hydrodynamic instabilities are known 
to occur when the Reynolds number becomes large [231 178] . In the case of Oseen vortices, we 
show in Section [4.2l that a rapid rotation (i.e., a large circulation number) has a stabilizing effect 
on the vortex : the size of the local basin of attraction increases, and the non-radially symmetric 
perturbations have a faster decay as t — > oo. These empirically known facts can be rigorously 
established, although optimal spectral and pseudospectral estimates are not available yet. 

In the final section of these notes, we consider the particular situation where the initial 
vorticity is a superposition of iV Dirac masses (or point vortices). The corresponding solution of 
the two-dimensional vorticity equation is called the viscous N -vortex solution, and the goal of 
Section [5] is to investigate its behavior in the vanishing viscosity limit. Our main result, which 
is stated in Section [5. 1[ asserts that the viscous N- vortex solution is nicely approximated by a 
superposition of Oseen vortices whose centers evolve according to the Helmholtz-Kirchhoff point 
vortex dynamics |57[ [62] . This approximation is accurate as long as the distance between the 
vortex centers remains much larger than the typical size of the vortex cores, which increases 



through diffusion. In Section 15.21 we decompose the viscous iV-vortex solution into a sum of 
Gaussian vortex patches, and we introduce appropriate self-similar variables which allow us to 
formulate a stronger version of our result, taking into account the deformation of the vortices 
due to mutual interactions. The proof involves many technical issues which cannot be addressed 
here, but we sketch the main arguments in Section 15.31 and refer the interested reader to [28] 
for more details. In particular, we show in Section 15.31 how to systematically construct an 
asymptotic expansion of the viscous iV-vortex solution, and we briefly indicate how the error 
terms can be controlled once a sufficiently accurate approximation is obtained. 

The content of the present notes is strongly biased toward the scientific interests of the au- 
thor, and does not provide a comprehensive survey of all important questions in two-dimensional 
fluid dynamics. We chose to focus on self-similar vortices, but other types of flows such as vortex 
patches or vortex sheets also lead to interesting and challenging problems, especially in the van- 
ishing viscosity limit. Also, we should keep in mind that all real fluids are contained in domains 
with boundaries, so a comprehensive discussion of two-dimensional fluid mechanics should cer- 
tainly include a description of the flow near the boundary, a question that is totally eluded here. 
On the other hand, we do not claim for much originality in these notes : all results collected 
here have already been published elsewhere, although they were never presented together in a 
unified way. In particular, most of the results of Sections [3] and H] were obtained in collaboration 
with C.E. Wayne [301 131 j . and the content of Section [5] is entirely borrowed from [28]. The 
material presented in Section [2] is rather standard, although our proof of Theorem [23] is perhaps 
not explicitly contained in the existing literature. The uniqueness part in Theorem 12.81 which 
is briefly discussed at the end of Section 13.31 was obtained in collaboration with I. Gallagher 
[25] . Finally, in Section [4.21 our approach to study the properties of the linearized operator at 
Oseen's vortex in the large circulation regime was developed in a collaboration with I. Gallagher 
and F. Nier [26] • 

2 The Cauchy Problem for the 2D Vorticity Equation 

We consider the two-dimensional incompressible Navier-Stokes equations : 

dtu(x,t) + (u(x, t) ■ V)u(x,t) = uAu(x,t) Vp(x,t) , 

P (2.1) 

div u(x,t) = , 

where x S M 2 is the space variable and t > is the time variable. The unknown functions 
are the velocity field u(x,t) = (ui(x,t),U2(x,t)) G M 2 , which represents the speed of a fluid 
particle at point x and time t, and the pressure field p(x, t) S R. Eq. (|2.ip contains two physical 
parameters, the kinematic viscosity v > and the fluid density p > 0, which are both assumed 
to be constant. 

Equivalently, the motion of a planar fluid can be described by the vorticity field : 

u>(x,t) = diu 2 (x,t) - d2Ui(x,t) , 

which represents the angular rotation of the fluid particles. If we take the two-dimensional curl 
of the first equation in (|2.ip . we obtain the evolution equation 



d t LL>(x,t) + u(x,t) -Vuj(x,t) = uAco(x,t) , (2.2) 

which is the starting point of our analysis. Note that u ■ Vw = div(uw) because diva = 0. In 
the case of a perfect fluid (v = 0), Eq. (J2.2J) simply means that the vorticity is advected by the 



velocity field u(x,t) like a material particle. For real fluids (v > 0), the vorticity also diffuses at 
a rate given by the kinematic viscosity. 

The vorticity equation (|2.2p is definitely simpler than the original Navier-Stokes system (I2.ip . 
but it still contains the velocity field u(x,t). To make (|2.2p independent of (|2.ip . it is possible 
to express u in terms of oj by solving the elliptic system 

diui + d 2 U2 = , d\U2 - d 2 ui = ui . (2.3) 

If oj decays sufficiently fast at infinity (see Lemma 12.11 below), the solution is given by the 
two-dimensional Biot-Savart formula : 



«CM) - ttz / .1 .. 2 "(i/.<)dy. ( 2 - 4 ) 



J_ f ( x -yY 

27r J R 2 \x - y\ 

where x 1 - = (— X2,x±) and \x\ 2 = xf + x 2 if x = {x\,X2) G K 2 . The vorticity equation (J2.2J) . 
supplemented with the Biot-Savart law (|2.4() . is now formally equivalent to the Navier-Stokes 
equations (|2.ip . Once a solution to (|2.2p . (|2.4j) is found, the pressure field can be recovered up 
to an additive constant by solving, for each t > 0, the Poisson equation 

— Ap = pdiv((u ■ V)u) , p{x,t) — > as \x\ — > oo . 

The following standard result shows that the Biot-Savart formula (12. 4p is well defined if the 
vorticity u lies in L P (M. 2 ) for some p G (1, 2). 

Lemma 2.1 Assume that oj G L p (M. 2 ) for some p G (1, 2). TTien i/ie velocity field u defined (for 
almost every x G M?) by the Biot-Savart formula (|2.4p satisfies : 
i) u G L 9 (IR 2 ), where q G (2,oo) is suc/i i/iai ^ = ^ — \- Moreover 

||«|U9(R2) < C|M| LP(R 2) . (2.5) 

mJ Vn G L P (R 2 ), and (H3D /joZds in L P (R 2 ). Moreover 

II v, "IIlp(r 2 ) < C|I w IIlp(ir 2 ) • ( 2 -6) 

Proof. According to f|2.4|) we have u = K * oj, where 

K( - x ) = T- f]2 > x G M 2 \ {0} , 

27T |xp 

is the Biot-Savart kernel. Since K belongs to the weak L 2 space L 2 '°°(IR 2 ), the Hardy-Littlewood- 
Sobolev (or weak Young) inequality [lH Section 4.3] shows that 

IMIl<2(r2) = ll-RT * c^||z,9(r 2 ) < 11-^ Hz, 2 .°°(r 2 )IMIlp 



ifl<p<2<oo and - = ^. This proves (|2.5p . Moreover Vn = X7K * oj, where \7K is an 

integral kernel of Calderon-Zygmund type which is associated to a bounded Fourier multiplier. 
It follows [721 Section 1.5] that K defines a bounded linear operator in L P (R 2 ) for p G (1, 2), as 
asserted in (12. 61). □ 



Remark 2.2 If oj G lJ-(WL 2 ), the velocity field u = K*oj belongs to L 2 '°°(]R 2 ), butu £ L 2 (R 2 ) in 
general, see \3J$ - If oj G L 2 (IR 2 ) ; one can solve (J2.3D directly by a simple calculation in Fourier 
space, which shows that u belongs to the homogeneous Sobolev space H^-iR 2 ), but u $ L°°(R 2 ) 



in general. Finally, if uj £ L P (M?) for some p £ (2,oo) ; the solution of (|2,3p is no longer given 
by the Biot-Savart law ()2.4p , oui 6?/ a modified formula of the form 

In t/iat case, u(x) may grow Zifce |x| *> as \x\ — > oo, and the solution of (J2.3J) is on/y defined up 
to an additive constant. 

2.1 General Properties of the Vorticity Equation 

In this section, we list a few important properties of the vorticity equation (|2.2p . including 
conservation laws, Lyapunov functions, and scaling invariance. The presentation is formal in the 
sense that we assume here that we are given a solution of (|2,2p with the appropriate smoothness 
and integrability properties. However, all calculations can be justified once suitable existence 
theorems have been established [52J. 

2.1.1 Conservations Laws 

Let uj(x,t) be a solution of (|2,2|) with initial data uo(x) = u(x,0). If ujq £ L 1 (1R 2 ), then 
x i— > uj(x, t) is integrable for all t > and the total circulation is conserved : 



oj(x,t)dx = / u>o(x)dx . (2-7) 

JR 2 J R 2 

If moreover \x\ujq £ -L 1 (M 2 ), then the first order moments of uj{x, t) are also conserved : 

XiUi(x,t)dx = / XiU)o(x)dx , i = 1,2 . (2-8) 



Finally, if in addition |x| 2 u;o £ L 1 (M 2 ), then the symmetric second order moment of io(x,t) 
satisfies : 

/ \x\ 2 uj(x,t)dx = / \x\ 2 u>o(x)dx + Aut I u)o{x)dx . (2.9) 

jr 2 Jm 2 Jm 2 

Thus the symmetric moment is conserved for perfect fluids (y = 0), and for viscous fluids (u > 0) 
if the total circulation vanishes. These elementary properties are easily established by direct 
calculations, using (|2.2p and the Biot-Savart law u = K * uj. In particular, to prove (j2.8|) . one 
uses the fact that 

/ Ui{x)ui(x) dx = K{(x — y)uj{y)uj{x) dx dy = , i = 1, 2 , 

7R 2 JR 2 xR 2 

because the Biot-Savart kernel K is odd. Similarly, the proof of (|2.9p relies on the identity 

/ (x ■ u(x))oo(x) dx = - I \{x — y) ■ K(x — y)\uj(y)uj{x)dxdy = 0, (2-10) 

JR 2 2 7 R 2 xK 2\ / 

which holds because K is odd and x ■ K(x) = for all x £ M. 2 . 
Remark 2.3 Sy Green's theorem, the total circulation satisfies 

7 := / ojdx = lim ® (^i(x) dxi + U2(x) dx2) , 



«->o°J|a:|=fl 



hence 7 represents the circulation of the velocity field at infinity. This quantity is conserved 
due to the structure of Eq. (12, 2p . In contrast, the conservation of the first- order moments and 
the simple evolution law for the symmetric second order moment are related to translation and 
rotation invariance, and would not hold if the vorticity equation was considered in a nontrivial 
domain Q C M 2 . 

2.1.2 Lyapunov Functions 

If ujq G L P (R 2 ) for some p G [l,oo], the solution of (|2.2j) with initial data ujq lies in L P (R 2 ) for 
all times and 

\\^(-,t)\\ L p (m.2) < ||o;o[|iP(]R2) , for all t > . (2.11) 

This again follows from the fact that the advection field u in (I2.2p is divergence-free. If the 
vorticity equation is considered in a domain Hcl 2 , and if one assumes as usual that u = on 
dQ (no-slip boundary condition), then vorticity is created in the boundary layer near d£l and 
(J2TTJI does not hold. 

If no = K * ujq G L 2 (M. 2 ), the solution u(x, t) of the Navier-Stokes equation (I2.ip with initial 
data no satisfies u(-,t) G L 2 (E 2 ) for all times, and the kinetic energy 

E(t) = - I |n(x,t)| 2 dx , 



is a nonincreasing function of time : 

E'(t) = -v / |Vn(x,t)| 2 d:r < . 
Jr 2 

The kinetic energy is the most famous Lyapunov function in fluid mechanics, and often the 

only one that is available. In particular, it plays a crucial role in the construction of global 

weak solutions of the three-dimensional Navier-Stokes equations |47| . Somewhat paradoxically, 

in unbounded two-dimensional domains, this quantity is not always useful because it is finite 

only if the total circulation vanishes. Indeed, we have the following elementary result : 

Lemma 2.4 If u G L 2 (R 2 ) 2 satisfies u = d\U2 — d^U\ G L 1 ^ 2 ), then J R2 ujdx = 0. 

Proof. By assumption the Fourier transform ui(k) = ikiii,2(k) — ik2Ui(k) is a continuous function 
of k G M 2 , hence 



/ udx = o)(0) = lim — 7T / Cj(k)Ak 
Jr2 e->0 7re z J\u <e 



i\k\ 
But 

—J \Cj{k)\Ak < — 2 \ \ k \\u(k)\dk <-}==( [ \u(k)\ 2 dkV /2 , 
ne J\k\<e nt J\k\<e V2vr \J\k\<e ' 

thus taking the limit e — > we obtain J R2 oj dx = 0. □ 

As a substitute for the kinetic energy, one can consider the pseudo-energy 

£d(t) = ~r / lo §i ?u(x,t)u(y,t)dxdy , 

47r JR 2 xR 2 \x-y\ 

where d > is an arbitrary length scale. This quantity is well defined under mild integrability 
assumptions on w, and can be finite even if the total circulation is nonzero. Morever, £d(t) is a 
Lyapunov function : 

S' d (t) = -v I uj(x,t) 2 dx < , 



and it is straightforward to verify that £d(t) = E(i) if the velocity field associated to co(x,t) is 
square integrable. In this particular situation, Edit) does not depend on the parameter d. 



2.1.3 Scaling Invar iance 

In the whole plane R 2 , the Navier-Stokes equations (12. lh and the vorticity equation (12. 2h are 
invariant under the rescaling 

u(x,t) ^Xu(Xx,X 2 t) , uj(x,t) ^ X 2 uj(Xx,X 2 t) , (2.12) 

for any A > 0. To solve the Cauchy problem, it is rather natural to use space-time norms which 
are invariant under the transformation (12.121) . The corresponding function spaces are called 
critical spaces (with respect to the scaling). A typical example is the energy space 

u£X = C$(R + ,L 2 (R 2 )) , \\u\\ x = sup|K,t)|| L2 , (2.13) 

where C^ (R+,£ 2 (R 2 )) is the space of all bounded and continuous maps from R + = [0, oo) into 
L 2 (R 2 ), and u(-,t) denotes the map x t-t u(x, t). Alternatively to include solutions with nonzero 
total circulation, one can assume that the vorticity is integrable and use the space 

ueY = C b °(R+,L 1 (R 2 )) , ||w||y = sup||w(-,t)|| L i . (2.14) 

t>0 

A well-known "metatheorem" in nonlinear PDE's asserts that, for a nonlinear partial differ- 
ential equation with scaling invariance, the critical spaces are the largest ones (in terms of local 
regularity of the solutions) in which we can hope that the Cauchy problem is locally well-posed. 
In fact, there may be other obstructions to local well-posedness, but according to the above 
claim it is certainly not reasonable to try to solve the Cauchy problem in function spaces that 
are "strictly larger" than the critical ones. 

As far as the Navier-Stokes equation is concerned, we know since the work of Leray [46] 
that Eq. (|2.ip is globally well-posed in the energy space f|2.13|) . and a similar result holds for 
the vorticity equation (J2.2J) in the critical space (J2.14J) . see [31 [6]. In the rest of this section, we 
discuss the Cauchy problem for (|2.2p in the space of finite measures, which is a natural extension 
of L 1 (M 2 ). As we shall see in Section El this generalization is essential to study the interaction 
of vortices in the vanishing viscosity limit. 

2.2 The Space of Finite Measures 

Let J\A(M. 2 ) denote the space of all real- valued Radon measures on R 2 , equipped with the total 
variation norm : 



|/u|| tv = sup<^ / <f)dn 



eO,(K 2 ), \\<f>\\ LOB <l} . (2.15) 

Here Co(R 2 ) is the space of all continuous functions <\> : R 2 — > R such that <f>(x) — > as \x\ — > oo. 
It is well known that A^(R 2 ) equipped with the total variation norm is a Banach space, which 
contains L X (R 2 ) as a closed subspace (if we identify any absolutely continuous measure with its 
density function). In particular, if u G L X (R 2 ), then ||w||tv = IMIz, 1 - The total variation norm is 
scale invariant in the sense that ||/x,\||tv = IIHItv f° r an ^ > 0, where fi\ is the rescaled measure 
defined by 

/ 0d/i A = / (piX^^dfi , for all <p £ C (R 2 ) . 

As is clear from the definition (|2.15p . the space A4(R 2 ) is the tolopogical dual of Co(R 2 ). By 
the Banach- Alaoglu theorem, it follows that the unit ball in A^(R 2 ) is a (sequentially) compact 
set for the weak convergence defined by : 

[i n v /i if / (l)dn n > / 4>dfjL for all <\> G C (R 2 ) . 



Any n £ Ai(M. 2 ) can be decomposed in a unique way as \x = f/, ac + fJ, s , where fi ac £ Ai(M. 2 ) is 
absolutely continuous with respect to Lebesgue's measure on R 2 , and fj, s is singular with respect 
to Lebesgue's measure [67] . By the Radon-Nikodym theorem, there exists a unique ui £ L X (R 2 ) 
such that dfi ac = udx. Furthermore, the singular part /u s can be decomposed as fi s = /i sc + /i pp , 
where \i vv is the restriction of [i to the atomic set £ = {x £ R 2 \ fi({x}) ^ 0}. Since fi is a finite 
measure, this set is at most countable, so that 

oo 
VPP = ^2li 6 Xi > ( 2 - 16 ) 

for some "fi £ R and some x« £l 2 . By construction, the singularly continuous part /jL sc has no 
atoms, yet is concentrated on a set of zero Lebesgue measure. 

To summarize, we have [i = [i ac + [i sc + [i pp , with [i ac _L [i sc _L /i pp . This notation means that 
the absolutely continuous, singularly continuous, and pure point parts of \i are concentrated on 
pairwise disjoint sets. It follows in particular that 

oo 

IIHItv = H/^acHtv + ||Msc||tv + H^pplltv = II^IIl 1 + H/^sclltv + 2-/ l7i| • 

i=\ 

If we now suppose that the measure fi £ Ai(M?) is the vorticity of a two-dimensional flow, 
there is a nice correspondence between the abstract decomposition result presented above and 
the following standard catalogue of nonsmooth flows : 

• Vortex Patches. A vortex patch can be defined as a flow for which the vorticity \x is a 
piecewise smooth function. For instance, /i can be the characteristic function of a smooth 
bounded domain Q, C M. 2 . In that case \x = fi ac £ L^IR 2 ) n L°°(1R 2 ), and the velocity field 
u = K * [x is Lipschitz continuous if the boundary dVl is sufficiently smooth. 

• Vortex Sheets. In contrast, vortex sheets are characterized by discontinuities of the velocity 
field. As a typical example, we mention the case where the vorticity [i is the Euclidean 
measure supported by a smooth closed curve 7 C M 2 . In that case \x = fi sc , and the 
tangential component of the velocity field u = K * /j, is discontinuous on the curve 7. 

• Point Vortices. This is the case where the vorticity /i is a collection of Dirac masses. Here 
/j, = fipp, and the velocity field u = K * /x is unbounded near the center of each vortex. 

2.3 The Cauchy Problem in M(R 2 ) 

The aim of this section is to study the Cauchy problem for the two-dimensional vorticity equation 
(J2.2I) in the space of finite measures. We start with a preliminary result that is relatively easy 
to prove. 



Theorem 2.5 [HJ | 

There exists a universal constant Cq > such that, if the initial vorticity fj, £ .M(R 2 ) satisfies 
W^ppWtv < Cqv , then the vorticity equation (|2.2|) has a unique global solution 

uj£C°((0,oo),L 1 (R 2 )nL co (M 2 )) 

such that \\uj(-,t)\\ L i < ||/x|| tv for all t > 0, and u(-, t) —^/i as t —$■ 0. 



Existence of global solutions to (|2.2p with initial data in A4(M?) was first proved by G.- 
H. Cottet [E], by Y. Giga, T. Miyakawa and H. Osada [M], and by T. Kato [43J. In these 
works, the strategy is : i) to approximate the initial measure by a smooth vorticity distribution; 
ii) to construct a global solution, with good a priori estimates, by applying classical existence 
results for smooth initial data; hi) to obtain a solution of the original problem by extracting 
an appropriate subsequence as the regularization parameter converges to zero. This yields the 
existence part of Theorem l2.5l without any smallness condition on the initial data. The restriction 
||/%>||tv < Co v arises when one tries to prove uniqueness of the solution by a standard application 
of GronwalPs lemma. As we shall see later, this smallness condition is purely technical and can 
be completely relaxed, see Theorem 12.81 

In the rest of this section, we present a more direct proof of Theorem[23J which is nevertheless 
inspired by [34] . Without loss of generality, we assume from now on that the kinematic viscosity 
is equal to 1. Given initial data /i G .M(M 2 ), we consider the integral equation associated with 
(1231) : 

u(t) = e tA fi - J div(V*- s ) A «(a)w(s)) ds , t > , (2.17) 

where ui(t) = oj(-,t), u(t) = u(-,t), and e denotes the heat semigroup defined by 

(e tA fi)(x) = -i- / e -l^l 2 /(4t) dfi y , t > , xeR 2 . (2.18) 

4vrt J R 2 

Our goal is to solve the integral equation (12.170 by a fixed point argument in an appropriate 
function space. This is a classical method which, in the context of the Navier-Stokes equations, 
goes back to the work of Kato and Fujita [21]. Solutions of the integral equation (|2.17p are 
often called mild solutions of the vorticity equation (|2.2p . It is rather straightforward to show 
that, if u G C°((0,T),L p (R 2 )) is a mild solution of (f2~21) with p > 4/3, then u(x,t) is smooth 
and satisfies (|2,2p in the classical sense. Thus there is no loss of generality in considering ([2.17P 
instead of (|231) . 

2.3.1 Heat Kernel Estimates 

The following lemma plays a crucial role in the proof of Theorem 12.51 : 

Lemma 2.6 Let fx G M(M 2 ). 

a) For 1 < p < oo and t > 0, we have 

||e tA /u|| LP < — -IIHItv , \\Ve tA fi\\ LP < - TZI M\tv ■ (2.19) 

(4-Kt) p t 2 p 

b) For 1 < p < oo, we have 

L p (/j,) := limsup(47rt) ~p\\e tA fi\\LP < \\fipp\\tv ■ (2.20) 

t-¥0 

Proof. Estimates (J2.19P are easily obtained from (|2.18p if p = 1 or p = oo, and the general 
case follows by interpolation. Thus we concentrate on estimate (|2.20p . which we establish using 
a simplified version of the argument given by Giga, Miyakawa, and Osada [34J. First, if we 
decompose p, = fj, ac + fi sc + /i pp , we have 

Lp{p) < L p (fi ac ) + L p (fj, sc ) + L p (fip P ) < Lp(n ac ) + Lp(n sc ) + ll/Upplltv ■ 
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Thus, to prove (|2.20p . it is sufficient to show that L p (/j,) = for all \x G Ai(M. 2 ) such that p pp = 0. 
Furthermore, since 

L p (jj) < LyipflPLM 1 - 1 '* < llMll^LooC/i) 1 - 1 /? , 

it is sufficient to consider the case where p = oo. 

Assume now that p G Ai(M?) satisfies p pp = 0, and fix e > 0. We claim that there exists 
5 > such that 

sup \p\ (B(x,S)) < e , where B(x,S) = {y G M 2 \ \y - x\ < 5} . (2.21) 

Here and in what follows, we denote by \p\ G A4(M. 2 ) the total variation of the measure p G 
A^M 2 ). If p + ,p,~ are the positive and negative variations of p, then p = p + — p~ (Jordan's 
decomposition) and \p\ = p + + p~ [67]. To prove (|2.2ip . assume on the contrary that there exist 
8 n > and x n G M 2 such that 5 n — > as n — > oo and \p\(B(x n ,6 n )) > e for all n G N. Since \p\ is 
a finite measure, the sequence (x n ) is necessarily bounded, thus after extracting a subsequence 
we can assume that x n — > x as n — > oo, for some i£l 2 . As a consequence, for any 5 > 0, we 
have B(x,5) D B(x n ,S n ) for all sufficiently large n, so that \p\(B(x,5)) > \p\(B(x n ,5 n )) > e. 
Taking the limit 5 — > 0, we obtain |//({x})| > e, which contradicts the assumption that p pp = 0. 
This proves (12311 . 

Now, for any t > 0, we can take x(t) G M. 2 such that |(e tA /x)(x(£))| = ||e iA ^||£°o. Thus, using 
flZHSD and (123TD . we obtain 

47rt||e /i||Loo < / e 4* d|/% + / e « d\p\ y . 

JB(x(t),5) JB{x{t),&Y 

The first term on the right-hand side is bounded by e for all t > 0, and the second one vanishes 
as t — > 0. Thus 

Loo(p) = limsup (4irt) ||e* mIU 00 < e • 

Since e > was arbitrary, we conclude that L OQ (p) = whenever // pp = 0, which is the desired 
result. □ 

Remark 2.7 Estimate (|2.20p was strengthened by Koto JJW m i/te following way : 

1 oo 

lim(47rt) 1 -^||e 4A /i|| L P = p~ 1,p \\{li}T=i\W < £ N , (2-22) 

fr->U i=l 

where ji , x% are as in (J2.16P . In 6oi/i (|2.20p , (|2.22j) ; the assumption that p > 1 is crucial : if 
p G M(M. 2 ) is a positive measure, then \\e t p\\ii = \\p\\t v for all t > 0, hence L\(p) = ||/u|| tv . 

2.3.2 Fixed Point Argument 

To prove the existence of solutions to the integral equation (|2.17p . we fix T > and introduce 
the function space 

X T = {w G C°((0,T],L 4 / 3 (IR 2 )) I \\u\\ Xt < oo} , (2.23) 

equipped with the norm 

||w||x T = sup t 1/4 ||u;(i)|| L 4/3 • 

0<t<T 
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Given \i £ Ai(M. 2 ), we denote u>o(t) = e [i for all t > 0. In view of Lemma [2.6l we have wo £ -Xt 
and there exist positive constants C\,C2 (independent of T) such that 

a) 1 1 wo \\x T < Ci\\fi\\ tv (for any T > 0); 

b) ||wo||x T < C2 1 1 A*pp 1 1 tv + e (if T > is sufficiently small, depending on fj,,e). 
Here e > is an arbitrarily small positive constant. 

On the other hand, given ui G Xt, we define Fuj : (0, T] — > L 4 / 3 (M 2 ) by 

(Foj)(t) = f div(e (t ~ s)A u(s)u(s)\ ds , < t < T , (2.24) 

where u(s) = K * u>(s) is the velocity field obtained via the Biot-Savart law (|2.4h . Using the 
second estimate in (|2.19p . Holder's inequality, and estimate (|2.5p . we find 

t^\\(FLo)(t)\\ Liri < &<■ f —^ \\u(s)u(s)\\ L i ds 

JO (t — SJi 

< t l/A f —^\\u{s)\\ L 4u{s)\\ Li/3 ds 

Jo (t-s)i 

<t l/A f— ^-j\\co(s)\\ 2 L4/3 ds (2.25) 

JO (t-S)i 

< Cll^llLi 174 / , C 3 i ds < C\\ U \\ 2 Xt , 



'0 (t — s)iS2 

for all t £ (0, T]. This shows that ||Fa;||x T < C3 IMI 2 ^ for some C3 > (independent of T), and 



it is not difficult to verify that (Fu)(t) depends continuously on t in the topology of L 4 ' 3 (R 2 ), 
so that Fui E Xt- Similarly, one can show that 

\\Foj - FuJ\\x T < C 3 (\\u\\x T + ||w||x T )||w - u\\x T , (2-26) 

for all u,u> € Xt- 

Now, fix R > such that 2C3R < 1, and consider the closed ball 

B = {ujEX t \ \\uj\\ Xt <R}CX T - 

As a consequence of the estimates above, if [|o;o||x T < -R/2, the map u; 1— >• wo — -Fu; is a strict 
contraction in B, hence has a unique fixed point in B. By construction, this fixed point uj is 
a solution of the integral equation (|2. 1T[) on (0, T]. The condition ||wo||x T < -R/2 is fulfilled if 
either 

1) 2Ci||//||tv < -R and T > is arbitrary, or 

2) 2C2||/i pp ||tv < R and T > is small enough, depending on /i. 

In other words, the proof above shows that equation (J2.17P is globally well-posed for small 
initial data \x £ .M(K 2 ), and locally well-posed for large initial data with small atomic part /i pp . 
However, if 2C2||/%>||tv > R, it is not possible to meet the condition [|wo||x T < -R/2 by an 
appropriate choice of T, so the argument breaks down and does not give any information on the 
existence of solutions to (|2.17p . 
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2.3.3 End of the Proof of Theorem [231 

The solution w G Xy of (12.17P constructed by the fixed point argument automatically satisfies 
w G C°((0,ri,L 1 (R 2 )). Indeed, arguing as in l|235j> . we find 

||w(t)-w (*)[Ui < / — ^-rl|o;( S )||2 4/ 3d S < C[M& T , (2.27) 

Jo (t-s)2 

for t G (0, T], where ^o(^) = e /i. Since H^oWIIl 1 — ll/^lltv by (12.19p . we conclude that 
ui(t) G L 1 (M 2 ) for all £ G (0,T], and the continuity with respect to time is again easy to verify. 
In fact, one can even show that 

limllwft) -w (t)||ii = . (2.28) 

t->o 

Since uJo(t) — *■ /x as t — > 0, this implies that w(t) — ^ \i as i — )• 0. Moreover, as ||uj(i) ||_t,i is a 
nonincreasing function of t, we deduce from (|2.28j) that ||u;(£)|| L i < ||/x|| tv for all t G (0,T]. 
To prove (|2.28|) . we denote 

5 := hmsup t 1 ' \\oj(t) — u>o(t)\\ L 4/3 = limsup \\u — uq\\x t ■ 
t^o r->o 

Since ui — ujq = (Fojq — FoS) — Fojq and ||w||x T + ll w o||x T < 2i?, it follows from (|2.26f) that 
S< (2C 3 i?)<5 + 4 /3 ( ) u), where 

^p(aO = hmsupt p ||(Fo;o)w||lp , 1 < p < oo . 

t-»-o 

Now, a direct calculation, which is postponed to Lemma 12.91 below, shows that £ p (p) = for 
any \i G Ai(M. 2 ) and any p G [l,oo]. This implies that 5 = because 2C3R < 1. Using again 
the identity uj — uiq = {Fuiq — Fuj) — Fooq and arguing as in (|2.27|) . we conclude that 

limsup \\u{t) -ua{t)\\iA < CR6 + £i([j,) = 0. 

To finish the proof of Theorem 12. 5| it remains to show that all solutions of (|2.2p are global 
in time. We first observe that, if \i G M(M?) n L 4 / 3 (R 2 ), then \\e tA fi\\ L 4/ s < \\fJ-\\ L 4/z for all 
t > 0, hence ||wo||x T ^ ^ 1 II^IIl 4 /3- Since the local existence time T > is determined by 
the condition ||wo||x T < -R/2, we see that an upper bound on the L 4 ' 3 -norm of the initial data 
provides a lower bound on the local existence time T. Assume now that fi G .M(R 2 ) satisfies 
2C 2 ||//pp||tv < R, and let U) G Xt be the local solution of (12. 2p constructed by the fixed point 
argument. If T\ > satisfies T x ||u;(T)|| L 4/3 < R/2, the same argument provides a solution cj 
of (|2.2p on the time interval [T, T + Ti] with Co(T) = U)(T), and arguing as in [3] one can verify 
that u G C°([T,T + T 1 ],L 1 (]R 2 )nL 4 / 3 (]R 2 )). Thus, if we glue together u and u, we obtain a local 
solution u G C°((0, T + Ti], L X (R 2 ) n L 4 / 3 (R 2 )) of fl22J). Iterating this procedure, and using the 
fact that ||w(t) 11^4/3 is a nonincreasing function of time, see (|2.1ip . we can construct a (unique) 
local solution on the time interval (0, T + kT\\ for any k G N. This means that the solution of 
(JZ2D is global if 2C 2 ||/x pp || tv < R. 

The proof of Theorem 12.51 is now complete except for two minor points. First, once we know 
that the solution u of (12. 2ft lies in L 4 ' 3 (R 2 ) for positive times, a standard bootstrap argument 
shows that u(t) G -L°°(R 2 ) for t > 0, as asserted in Theorem 12.51 More importantly, while 
the fixed point argument only shows that the solution is unique in a ball of the space Xt, one 
can in fact prove uniqueness in the whole space C 6 °((0, T], L 4 (R 2 )) n C°((0, T], L°°(R 2 )), for any 
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T > 0. Establishing this improved uniqueness property requires additional arguments which 
will be presented in Section T3.ll below. 

Theorem 12.51 shows that the vorticity equation (|2.2p is globally well-posed if the atomic part 
of the initial data satisfies ||// p J|tv < Cqu for some Cq > 0. As was already mentioned, this 
smallness condition can be completely relaxed, and the optimal result is : 

Theorem 2.8 [MI55] 

For any fi £ 7W(M 2 ), the vorticity equation (|2.2[) has a unique global solution 



u 6 C°((0, oo), L 1 (R 2 ) n L°°(R 2 )) 
such that \\u)(-,t)\\ L i < ||/x|| tv for all t > 0, and u>(-,t) — *• // as £ — >■ 0. 

The existence part in Theorem 12.81 can be established using a relatively standard approxi- 
mation procedure, see [54] , but proving uniqueness without any restriction on the initial data 
requires a careful treatment of the large Dirac masses in the initial measure. We refer to Section[3] 
below for a sketch of the uniqueness proof, and to [25] for full details. 

2.3.4 Short Time Self-Interaction 

To conclude this section, we state and prove an auxiliary result which was used in the proof of 
Theorem 12.51 

Lemma 2.9 Let fi G A4(R 2 ). For any t > 0, let a>o(t) = e tA fi and let uo(t) = K * w (t) be the 
corresponding velocity field, given by the Biot-Savart law (12. 4\\ . Then 

l p (n) ■= limsup i 1_ p||(Fw )(t)||iP = , 1 < p < oo , (2.29) 

where 

(Foj ){t) = I e {t ~ s)A u (s) ■ Vw (s) ds , t > . 

Jo 

Proof. We shall prove (|2.29[) for p £ [1,2). The other values of p can be treated in a similar 
way. Arguing as in (|2.25p . we find 

^"pIKFwoXtJHtP < t l -'v [ — ||wo(s)|| 2 4/3 ds , t > . 

Jo (t-s) 2 £ 

It follows that £ p (n) < C p (L 4 / 3 (/j,)) 2 , where C p is a constant depending only on p and L p (fi) is 
defined in (|2.20p . This already shows that £ p {ijl) = if fi pp = 0. More generally, the argument 
above implies that 

Mm + /u 2 ) - £p(jJti) - i P (p2)\ < 2C p L 4/3 ( Pl )L 4/3 (fi 2 ) , (2.30) 

for any fii, 1x2 € .M(R 2 ). Taking /ii = \i — fi pp and fi 2 = ^ PP , we conclude that t p {fj) = £ p (/j, pp ). 
This means that the non- atomic part of the measure /U does not contribute to the limit £ p (fJ>)- 
It is thus sufficient to prove (|2.29|) in the case where \i is purely atomic, and in view of (|2.3U|) 
we can even suppose that fi is a finite linear combination of Dirac masses : \x = ^2i=i a A* • m 
that case, 



N N 

uo(x,t) = ^2t G ( — 7T) ' M x ,t) = Yj 



Oii /-1 / X X 



V 



,t-\ y/i J' uv ' ; {rtVt V vt 
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where G and v G are defined in (13.21) below. Let 



H(x,t) = u (x,t).Vu;o(x,t) =J2^v G (^-).VG(^-^- 

Here the sum runs over indices i ^ j because the contributions of the self-interaction terms i = j 
are identically zero. Using this observation, it is rather straightforward to verify that 

limsupt 2 f \\H(t)\\up < oo , 1 < p < oo . 

i-»0 

As p < oo, it follows that 

t l " P \\(Fu )(t)\\ LP < t 1 ^ I C L^\\H{s)\\ L ,ds —> . 

Jo (t-sf'p *- >0 

This concludes the proof. D 

3 Self-Similar Variables and Long-Time Behavior 

Explicit examples of two-dimensional viscous flows are easily constructed if one assumes radial 
symmetry. Indeed, if the vorticity u(x, t) is radially symmetric, the velocity field u = K*uj given 
by the Biot-Savart law (|2.4p is azimuthal, in the sense that x ■ u(x, t) = 0. As a consequence, the 
nonlinearity u ■ Voj in (|2.2[) vanishes identically, and (|2.2p therefore reduces to the linear heat 
equation dtOJ = vAlj (which preserves radial symmetry). 

As an example, consider the particular case where the initial vorticity is a point vortex 
of circulation 7 G K located at the origin: p, = 7^0- The unique solution of (12. 2D given by 
Theorem 12.81 is the Lamb- 'seen vortex: 

^ t » = ^ G (^)- -<*•'> = tfAtz) • (3 ' 1) 



where the vorticity and velocity profiles are given by 

OK) = JL«-KIV\ „o K ) = i-| 5 (l-eH«4), {eR2 . (3 . 2) 

Note that J R2 G(£) d^ = 1, so that 7 = J R2 u(x,t) dx is the total circulation of the vortex. 
Oseen's vortex is thus a self-similar solution of (|2.2p with Gaussian vorticity profile G. The 
velocity profile v is azimuthal, vanishes at the origin, and satisfies \v (£)| ~ (27r|^|) _1 as 
|£| — > 00. In particular, v G ^ L 2 (M. 2 ), in agreement with Lemma 12.41 Oseen vortices are 
therefore infinite energy solutions of the two-dimensional Navier-Stokes equation. 

Oseen's vortex plays an important role in the dynamics of the vorticity equation (|2.2p . As 
was already observed, it can be considered as the "fundamental solution" of (J2.2H : i.e., the 
solution with a single Dirac mass as initial data. In addition, the following result shows that it 
describes the long-time asymptotics of all integrable solutions of (|2.2p : 



Theorem 3.1 [331 [H] 

For all initial data fi 6 Ai(R 2 ), the solution u:(x,t) of (|2.2p given by Theorem \2.8\ satisfies 



lim 

t— >oo 



uj(xA) G[—=) = , where 7 = / da . (3.3) 
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It follows immediately from (|3.3p that Oseen vortices are the only self-similar solutions of 
the Navier-Stokes equation in R 2 for which the vorticity profile is integrable. As an aside, 
we mention that, following the approach of Cannone and Planchon [9|, one can construct an 
infinite-dimensional family of small self-similar solutions of (|2.ip for which u G L 2,00 (IR 2 ) but 
co $ L 1 (M 2 ). Assuming, as in Theorem 13.11 that the initial measure is finite allows to eliminate 
all these "artificial" solutions, which are produced by fat tails in the initial data and not by 
the intrinsinc dynamics of Eq. (|2.2p . Also, Theorem 13.11 strongly suggests that Oseen vortices 
are stable solutions of f)2.2f) for all values of the circulation 7. This important question will be 
discussed in Section H] below. 

It is worth mentioning that the convergence result (|3.3p is not constructive, and does not 
provide any estimate of the time needed for the solution of (|2.2p to approach Oseen's vortex. 
Under strong localization assumptions, such estimates were obtained in [29 , but these results 
are certainly not optimal. We also refer to [8] for an interesting attempt toward a more precise 
description of the intermediate asymptotics, using ideas from statistical mechanics. 

In the rest of this section, we give a sketch of the proof of Theorem 13.11 and of the uniqueness 
part of Theorem 12.81 We first recall classical estimates for solutions of convection-diffusion 
equations, which are due to Osada |63j . We then introduce self-similar variables, which allow 
to compactify the solutions of (12. 2p in L 1 (M 2 ), and to transform the self-similar Oseen vortices 
into equilibria. Finally, we use Lyapunov functions to determine the w-limit sets of solutions of 
(12. 2p . and we establish a "Liouville theorem" which implies both the convergence result (I3.3P 
and the uniqueness of the solution of (|2.2p with a single Dirac mass as initial vorticity. 



3.1 Estimates for Convection-Diffusion Equations 

Given v > 0, we consider a linear convection-diffusion equation of the form 

d t uo(x, t) + U(x, t) ■ Vw(i, t) = uAuj(x, t) , (3.4) 

where x G M 2 , t € (0, T), and U : M? x (0,T) — > M. 2 is a given divergence-free vector field. We 
assume that the curl Q = d 1 U 2 - d 2 Ui of the advection field U satisfies Q G C°((0,r),L 1 (R 2 )) 
and ||fi(-,£)|| L i( K 2) < Kqu for t G (0, T), where Kq is a positive constant. Then, in view of [63] . 
any solution of (|3.4j) can be represented as 

u(x,t) = / Tu(x,t;y,s)uj(y,s)dy , xeM 2 , <s <t <T , (3.5) 

where T^ is the fundamental solution of the convection-diffusion equation (|3.4j) . The following 
properties of Vfr will be useful: 

• There exist j3 > and K\ > (depending only on Kq) such that 

< m *,P,«,s) < -*Lj exp^J^Hl) , (3.6) 

for x, y G M 2 and < s < t < T, see [63]. In fact, a more precise result due to Carlen and Loss 
[TO] shows that, if {vt) l / 2 \\U{-, t)\\ Loa < K Q v for all t G (0,T), then d£5J holds for any G (0, 1), 
with a constant K\ depending on Kq and /3. It is also possible to establish a Gaussian lower 
bound on T^, see [63] , 

• There exists 7 G (0, 1) (depending only on Kq) and, for any 5 > 0, there exists K 2 > 
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(depending only on Kq and S) such that 

|r^(x,t;i/, a )-r^(x / ,t / ;i/,s / )| 

< k 2 (\x- x 'v + W-t')?' 2 + \y-vV + |v(«-a / )| 7/2 ) , (3.7) 

whenever v(t — s) > 5 and v(t' — s') > 5, see [63J. 
• For < s < t < T and x, y G R 2 , 

/ T^(x,t;y,s)dx = 1 , / T^(x,t;y,s)dy = 1 . (3.8) 

Jr 2 Jr 2 

For < s < r < t < T and x, y G M 2 , 

r^(x,t;y,s) = / ru(x,t;z,r)Tu(z,r;y,s)dz . (3.9) 

Jm 2 

Remark 3.2 If x,y G R 2 and £ > 0, it follows from (|3.7p t/iai i/ie function s i— >• r^(x, t;y, s) 
can fee continuously extended up to s = 0, and i/iat i/us extension (still denoted by Ty) satisfies 
properties (|3.6p to (|3.9p wra£a s = 0. 



As an application of these results, we give a new formulation of the uniqueness claims in 
Theorems [23] and ELl Assume that w G C ((0,T),L 1 (R 2 )nL oo (M 2 )) is a solution of ([12]) which 
is uniformly bounded in L 1 (IR 2 ) : there exists Kq > such that ||w(i)|| L i < Kqu for t G (0, T). 
Let 4> G C^°(M 2 ) be a test function. Using (|2.2p and integrating by parts, we easily find 

— / oj(x,t)4>(x)dx = / w(ar, t)u(x, t) ■ V(f)(x)dx + v / Lo(x,t)A(f>(x)dx . 

dt Jm.2 Jm.2 Jr2 

If we express u in terms of oj using the Biot-Savart law (|2.4p . we can rewrite the first integral 
on the right-hand side as 

1 f . ./_ *\. ./-.. ^ (a? ~ y) 1 - 

; xr P ~~ y\ 

see e.g. [20]. Since co is uniformly bounded in L 1 (M 2 ), we conclude that 



— / u(x, t)u(y, t) f y • (V^(x) - V0(y)) dx dy , 

47T ./iB2 v irp2 X-«r V / 



dt 



/ u(x,t)(f>(x)dx < Cv 2 V ||<9>|| L °o , te(p,T) , 

Jv 2 , , „ 



H=2 



for some constant C > depending only on Kq. This implies that oj(t) has a limit in D'(R 2 ) as 
£ — > 0, which we denote by /U. As ||w(t)||xi < Kqv, it follows that \i G A4(M 2 ) with ||//||tv < ^o^> 
and that w(i) — *• \i as i — >• 0. 

On the other hand, since w solves (|3.4|) with U(x,t) = u(x,t), and since the curl of the 
advection field is uniformly bounded in L 1 (M 2 ), we have the representation (13. 5h . where the 
fundamental solution F^(x,t;y, s) satisfies (|3.6p to (|3.9p . In particular, using Remark 13.21 we 
have for all xGK 2 and all t G (0, T), 



u(x,t) = / T^(x,t;y,0)uj(y,s)dy 

Jr 2 

+ / (r^(x, t; y, s) - T u u (x, t; y, 0))w(y, s) dy , < s < t 

Jr 2 
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In view of (|3.7p . the second integral on the right-hand side converges to zero as s — > 0. On the 
other hand, since y t-t r^(x, t; y, 0) is continuous and vanishes at infinity, and since U)(s) -^^as 
s — > 0, we can take the limit s — > in the first integral and we obtain the useful formula 

u(x,t) = / ri(x,t;y,0)diJ, y , i£l 2 , 0<£<T, (3.10) 

which shows that any solution u)(x,t) of (12.2p which is uniformly bounded in L 1 (IR 2 ) can be 
represented in terms of its trace \x at t = 0, using the corresponding fundamental solution r„. 

In particular, since T^(x,t;y,0) is positive and satisfies (I3.8p . it follows from (|3.10p that 
[|w(t)[|j^i < H/^Htv for all t E (0,T). Thus we can assume that Kq = ||//||tv without loss of 
generality. Moreover, using the upper bound (I3.6P we deduce from (I3.10P that 

sup (vtf'hwWWi* ^ Cll^lltv , 1 < V < oo , (3.11) 

0<i<T 

where C is a positive constant depending only on Kq. In fact, a more careful argument shows 
that the constant C > in (|3.1ip is independent of Kq, see [3l S21 [10] . Similarly, repeating the 
proof of Lemma 12.61 we find 

limsup(z4) p||cj(£)||.z> < C||/ip P ||tv , I < p < cc . (3-12) 

t->o 

Taking p = 4/3, we deduce from (|3.11|) . (|3.12p that the solution uj of (|2.2|) lies in the space 
Xt defined in (|2.23p . and that ||o;||x T is small if [|/i[|tv is small, or if ||/x pp ||tv and T are small. 
In both cases, u) necessarily coincides with the solution of (|2.2p constructed by the fixed point 
argument. 

3.2 Self-Similar Variables 

Let u e C°((0,oo),X 1 (M 2 ) nL°°(R 2 )) be a solution of (^j) . Given x E M 2 , i > 0, and T > 0, 
we introduce the self-similar variables 

* = "~ X ° , r = log(^) , (3.13) 

and we transform the vorticity and velocity fields as follows : 



1 ( x-xq , ft + t 

t + t \y/l/(t + t ) v T 



(3.14) 



. 1/ x-x ft + to 

u{x,t) = J— — - v log 



Then the rescaled vorticity w(£,t) satisfies the evolution equation 



d T w + v -V^w = A^w -\ — £ -V^w + w , (3.15) 



and the rescaled velocity v is again expressed in terms of w via the Biot-Savart law (|2.4p : 
v = i^ * u>. The initial data of (|2.2p and (|3. 15[) are related via 

w(£,t ) = t Q cj(x + ^y/l^,0) , £GM 2 , 
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where tq = log(to/T). Note that w(£,t) and v(£,t) are now dimensionless quantities, as are the 
new space and time variables £, r. If io(£, r) = aG(^) for some a£l, then by (13.141) u(x,t) is 
Oseen's vortex with circulation 7 = m/ located at point xo G M 2 and originating from time — to- 
Thus Oseen vortices are now equilibria of the rescaled system (|3.15j) . In what follows, we shall 
use the change of variables f)3. 14j) with xq = and to = T, so that tq = 0, but other choices of 
xq or to are sometimes more appropriate, see e.g. [29] . 

The main effect of the transformation (|3.14p is to replace the diffusion operator uA on the 
right-hand side of (J2.2I) by the more complicated operator 

C = A+-C-V + 1, (3.16) 

which appears in f|3. 15|) . As we shall see in Section^ the operator C has better spectral properties 
than the Laplacian, when acting on appropriate function spaces. Moreover, the semigroup e 
generated by C is asymptotically confining, a property that does not hold for the heat semigroup 
e ut . In fact, using (I2.18P and (|3.14p . it is straightforward to derive the explicit formula 

< e ^>«> = i^M I "K JS 7 a M /2 'V °'"» d " • < 3 - 17 > 

for all (£l 2 and all r > 0, where a(r) = 1 — e~ r . If wq G L 1 (R 2 ), then Lebesgue's dominated 
convergence theorem shows that (e T-C u>o)(£) converges to aG(£) as r — > 00, where a = L 2 wo d£ 
and G is given by (|3.2p . Thus the solution of the linear equation d T w = Cw is asymptotically 
confined like a Gaussian function provided the initial data are integrable. 

A similar property holds for the nonlinear equation (|3.15p . and allows to show that the 
trajectories of this system are compact. 

Lemma 3.3 For any wq G L 1 (M 2 ) ; the solution {w(t)} t >q of (|3.15p with initial data wq is 
relatively compact in L 1 (R 2 ). 

Proof. We first observe that the Cauchy problem for Eq. (|3.15p is globally well-posed in the 
space L 1 (M 2 ). This can be deduced, via the change of variables (|3.14p . from the corresponding 
statement for the original equation (|2.2|) . which is a particular case of Theorem 12.51 This 
argument also shows that the L 1 norm of the solution w(t) is nonincreasing with time, and 
that the circulation parameter a = f R2 w(£,,T~)dt; is a conserved quantity. Now, using the 
representation (I3.10p of the solution of (|2.2p in terms of the initial data, and the bound (I3.6P 
on the corresponding fundamental solution, we obtain the following estimate for the solution of 

KM| < 4r / exp(-/3 le 7 ? fip 2 )K(77)ld7 7 , (3.18) 

a {V 7r 2 v 4 a( T ) ' 

for ^6l 2 and r > 0, where again a(r) = 1 — e~ T . Here /3 G (0, 1) and K\ > depend only on 
IIwoIIl 1 - Using this bound, it is rather easy to show that 



sup/ K£,r)|d£- ►O, (3.19) 

T>1J\£\>R R^oo 

see |31l Lemma 2.5]. On the other hand, classical estimates for the derivatives of solutions of 
(|2.2p ensure that 

sup sup |Vm>(£j7")| < 00 , 

t>i§gr 2 
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see US 03]. This together with (l3TT9|) implies that 

sup sup / \w(Z, - n,r) -w(£,t) | d£ ► . (3.20) 

t>1 \ v \<8 JR 2 5 ^° 

By a well-known criterion of Riesz [66] . (|3.19|) and (|3,20|) together imply that the trajectory 
{w(r)} T >i is relatively compact in L 1 (M 2 ). Since we also know that w G C ([0, 1],L 1 (M 2 )), the 
conclusion follows. D 

If wq G L 1 (R 2 ), it follows from Lemma 13.31 that the solution w{r) of (J3.15P with initial data 
wq converges as r — > oo, in the L 1 topology, to the w-limit set 

n(w Q ) = f| Mr) I r > T} , 

T>0 

where the overline denotes here the closure in L 1 (R 2 ). This is a nonempty compact and connected 
subset of L 1 (M 2 ) which is positively and negatively invariant under the evolution defined by 
(J3.15J) . in the following sense: for any wo G £l(wo), there exists a complete trajectory w E 
C^L^R 2 )) of (f37T5]) such that w{t) G n(w ) for all r G 1 and w(0) = w . In the next 
section, we shall prove that Q(wq) consists of a single point {aG}, where a = L 2 ^o d^. 

3.3 Lyapunov Functions and Liouville's Theorem 

Lemma [3.3l shows that positive trajectories of (J3.15P in L 1 (M 2 ) are relatively compact. The same 
conclusion does not apply to negative trajectories in general. In fact, the following proposition, 
which is the main result of this section, implies that the equilibria associated to Oseen vortices 
are the only negative trajectories of A3. 15f) which are relatively compact in L 1 (M 2 ). 

Proposition 3.4 ("Liouville's theorem") If {u;(t)} tG r is a complete trajectory of (|3.15p which 
is relatively compact in L 1 ^ 2 ), then there exists a G R such that w(t) = aG for all r G R. 

Proof. Let A and Q denote the a-limit set and the w-limit set of the trajectory {w(t)} T £r in 

A = p| {w(t) \t <T} , n = f~] {w{t) I r > T} . 

T<0 T>0 

Both A, fl are nonempty compact sets in L 1 (R 2 ). Moreover, they attract w(t) in the sense that 
dist^i (w(t),A) — > as r — > — oo, while dist^i (w(t),Q) — > as r — > +oo. Following [31) . to 
characterize these limit sets more precisely, we use a pair of Lyapunov functions for Eq. (J3.15P . 

First step. Our first Lyapunov function is just the L 1 norm of the solution. Let <£ : L X (R 2 ) — > 
R+ be the continuous function defined by 



$H = / K£)|d£ , w£ L X (R 2 ) 

Jr 2 



and let 



In other words, a function w G L X (R 2 ) belongs to E if and only if u>(£) has almost everywhere a 
constant sign. We have already mentioned that $ is nonincreasing along trajectories of (J3.15D . 
as a consequence of (|2.1ip . Now, using the strong maximum principle, it is straightforward to 
verify that $ is in fact strictly decreasing, except along trajectories which lie in E : if {w(t)} t >q 
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is a solution of (|3.15j) such that $(w(t)) = $>(w(0)) for all r > 0, then w(0) G E (hence iy(r) € E 
for all r > 0). 

By LaSalle's invariance principle [71], the a-limit set of the trajectory {w(t)} t £r lies in 
the neutral set E, because A is positively invariant under the evolution defined by (J3.15P and 
entirely contained in a level set of the Lyapunov function <&. The same conclusion applies to 
the cj-limit set £1. In view of the definition of E, it follows that $(io) = |a| for all w £ A or £1, 
where a = / R2 w(£,t) d£ is the circulation of our solution (which is a conserved quantity). This 
in turn implies that &(w(t)) = \a\ for all rel, which is possible only if the whole trajectory 
{w(t)} tG m. lies in the neutral set E. Thus any relatively compact complete trajectory of (13.151) 
has necessarily a definite sign, which is the sign of the circulation parameter a. If a = 0, this 
means that w(£, r) = 0, in which case the proof is complete. If a < 0, we observe that the 
substitution w(£i,£2>t) *-* —w(£,2,£,i,t) leaves Eq. (|3. 15[) unchanged, but reverses the sign of 
the circulation. Thus we can assume henceforth that a > 0. In that case u>(£, r) > for all 
£ € R 2 and all r £ R, and proceeding as in Section [3.11 one can show that w(£,t) is bounded 
from above and from below by time-independent Gaussian functions. 

Second step. As a second Lyapunov function, we use the relative entropy of the vorticity 
distribution w(t) with respect to the Gaussian G [79]. If w : R 2 — > R+ is a (measurable) 
positive function with a Gaussian upper bound, we define 



H(w) 



w(?)log(|||df 



Then H is nonincreasing along trajectories of A3. 15[) . because 



— H(w(t)) 



w(£,t) 



Vlog 



G(0 



d£ < 



(3.21) 



(3.22) 



To prove f|3.22j) . we compute 

-^-H(w(t)) = [ (i + log^)d T wdt= [ (l + log^)(Cw-vVw)dC 

dr J R2 \ GJ J R 2\ GJ 

Using the identity Cw = div(GV(^)) and integrating by parts, we obtain 



/ (l + log^)(£u,)d£ = - / Vflog 



W 



W 



V(log^ 



G„/w\ 

w V (gH 



d£. 



On the other hand, using v ■ Vw = div(vw) and integrating by parts, we find 

/ ( 1 + log |J )(v ■ Vui) d£ = / (l + log(47™))(u-Vu;)df + / ^-(v-Vw)d£, 

JR2\ GJ J K 2 J R 2 4 



v ■ Vw d£ 



({,-v)wd£ = . 



Note that the last integral vanishes because v = K * w, see (|2.10p . This concludes the proof 
of inequality (13.220 . which shows in addition that H is strictly decreasing along trajectories 
of (|3.15p . except on the line of equilibria corresponding to Oseen vortices. More precisely, if 
{w(t)} t >o is a positive solution of (13.150 (with a Gaussian upper bound) such that H(w(r)) = 
H(w(0)) for all r > 0, then w(Q) £ T, where T = {/3G}/3>o is the half-line of positive equilibria. 
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We now return to the analysis of our complete trajectory {w(r)} Tg ]g of (|3.15j) . We know 
from the first step that w(t) G £ for all r£l, and without loss of generality we assume that 
the circulation parameter a = L 2 u>(£, r) d£ is positive, which implies that w(£, r) is positive 
and has a Gaussian upper bound. By LaSalle's invariance principle, both the ce-limit set A 
and the w-limit set f2 are contained in the line of equilibria T, because these sets are positively 
invariant and contained in a level set of the Lyapunov function H. The conservation of the total 
circulation then implies that A = fi = {aG}, hence H = alog(a) on A and £1. It follows that 
H(w(t)) = a log (a) for all t£|, which is possible only if w(r) G T for all r£l. But T is a 
line of equilibria, hence we must have w(t) = aG for all r G R, as stated in Proposition 13.41 □ 

As a first application of Proposition 13. 4} we show that Oseen vortices describe the long-time 
asymptotics of all solutions of (|2.2p in L 1 (M 2 ). 



Proof of Theorem SHI Assume that ui G C°((0,oo),L 1 (R 2 ) n L°°(R 2 )) is a solution of (1231) 
with initial data fj, G .M(R 2 ), in the sense of Theorem 1 2. 8\ and let w(£,, r) be the rescaled vorticity 
obtained from u;(x,t) via the change of variables (13.14p . with an arbitrary choice of the center 
x Q G R 2 and the initial time t = T > 0. Then w G C°((0,oo),L 1 (R 2 ) n L°°(R 2 )) is a solution 
of the rescaled vorticity equation A3. 15f) such that 



o 



/ w(^,r)d^ = — , where 7 = / w(x,t)dx = / d^ . (3.23) 

7R2 V J R 2 J R 2 



By Lemma 13.31 the trajectory {w(t)} t >i is relatively compact in L X (]R 2 ). Let ft denote the 
u;-limit set of this trajectory in L 1 (M 2 ). Then Q is a nonempty compact subset of L 1 (M 2 ) which 
is positively and negatively invariant under the evolution defined by (j3.15|) . In other words, Q 
is a collection of complete, relatively compact trajectories of f)3. 15|) . By Liouville's theorem, all 
such trajectories are equilibria of the form f3G, for some (3 G R. But the conservation of the 
total circulation implies that /3 = a, so we conclude that £1 = {aG}. This exactly means that 
w(t) — > aG in L 1 (M 2 ) as r — > 00. Returning to the original variables, we obtain (|3.3p . □ 

Remark 3.5 Under the assumptions of Theorem \3 '. 11 the above proof shows that 



lim 

t— >oo 



u(x,t)- , 7 t GJ * X ° =) =0, (3.24) 



for any xq G R 2 and any to > 0. This might be surprising at first sight, but there is no 
contradiction because, as is easily verified, two Oseen vortices with different values of xq or to 
converge to each other in the L 1 topology as t — >• 00. However, the choice of the vortex center 
and the initial time becomes important if one wants to determine the optimal decay rate in 
(|3.24p . This is possible if one assumes that the vorticity u(x,t) decays to zero sufficiently fast 
as \x\ — > oo ; see f31\ [2$/ and Section^ below. 

We next show that Proposition 13.41 implies the uniqueness of the solution of the vorticity 
equation (|2.2p with a point vortex as initial data. 



Proposition 3.6 [3D EZ] 

Let to G C°((0,T),L 1 (R 2 ) n L°°(R 2 )) be a solution of f$2?2$) which is bounded in L X (M 2 ) and 
satisfies u — v 7^0 as t — >• 0, for some 7 G R. Then 

u(x,t) = ^-G(-^=) , xGR 2 , 0<t<T. (3.25) 

vt V v vt ' 
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Proof. Applying (|3.10p in the particular case where // = 7^0, we obtain u(x, t) = 71^(2;, t; 0, 0), 
where T v u is the fundamental solution of (]3.4p with advection field u = K * uj given by the 
Biot-Savart law (12,41). In view of (13.61). we thus have 



|w(a?,<)| < ^1 exp(-p^-) , i£l 2 , 0<t<T. (3.26) 

for some K% > and /3 G (0, 1). We now apply the change of variables (|3.14p with xq = and 
to = 0. Since r = log(i/T), the rescaled vorticity w(£,t) is then defined for all £ E M? and all 
r < 0, so that w G C°((-oo,0),L 1 (R 2 ) n L°°(R 2 )) is a negative trajectory of (i3~T5|) . Moreover, 
the bound (1336]) implies that |w(£,r)| < {K l \^\/v)e-^ 2 / A , for all ?GK 2 and all r < 0. As in 
the proof of Lemma ET3] this bound implies that the trajectory {w(t)} t< o is relatively compact 
in L 1 ^ 2 ). Applying Proposition 13.41 we conclude that w(£,t) = aG(£) for some a G M, and 
the conservation of the total circulation implies that a = 7/z/. Thus w(£,t) = {^/v)G{£) for all 
£ G M 2 and all r < 0, and returning to the original variables we obtain (I3.25p . □ 

It is important to emphasize that there is no restriction on the circulation 7 in Proposition l3.61 
Thus the conclusion cannot be obtained by a standard application of Gronwall's lemma, which 
would require that the ratio \^)\/v be sufficiently small. Intuitively, Proposition 13.61 and Theo- 
rem [23] together imply that the solution of (|2.2|) with arbitrary initial data \i G Ai(M?) should 
be unique, as asserted in Theorem 12.81 Indeed, given any e > 0, we can decompose \x = \X\ + (J>2, 
where /j,% is a finite superposition of Dirac masses and || /-*2 lltv < e - If e is sufficiently small, we 
know from Theorem 12. 51 that fj2.2[) has a unique solution with initial data \x-i-, and Proposition l3.6l 
implies that each atom in \x\ also generates a unique solution of (12.2p . Although the the vorticity 
equation (|2.2p is nonlinear, this decomposition can be used to prove uniqueness if we can show, 
as in the particular case considered in Lemma 12.9] that the various components of the solution 
do not strongly interact for small times. These rough ideas can be turned into a rigorous proof, 
and we refer the interested reader to |25| for more details. 



4 Asymptotic Stability of Oseen Vortices 

As was mentioned in Section 13. 2| the rescaled vorticity equation 

d T w + v ■ V ? w = Cw , where £ = A+-f-V + l, (4.1) 

has a line of equilibria {aG} Q ^ which correspond, in the original variables, to the family of 
self-similar Oseen vortices (|3.ip . In fact, Proposition 13.41 implies that these are the only steady 
states of (|4.ip in the space L X (1R 2 ), and Theorem 13.11 even shows that the family {«G} ae s is 
globally attracting in the sense that any solution of (|4.ip in L 1 (M 2 ) converges to aG as r — > 00, 
where a = J R2 tf(^,r)d^ is the circulation parameter. This strongly suggests, but does not 
rigorously prove, that aG is a stable equilibrium of (|4.ip for any a£l. 

4.1 Stability with Respect to Localized Perturbations 

The aim of this section is to study the flow of (|4.ip in a neighborhood of Oseen's vortex aG, 
and to show that this steady state is indeed stable for any value of the circulation parameter 
aEl. Setting w = aG + w, v = av G + v, we obtain the perturbation equation 

d T w + v-Vw = (C- ak)w , (4.2) 

where 

Aw = v G ■ Vw + v ■ VG . (4.3) 
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Here v = K * w is the velocity field obtained from w via the Biot-Savart law (|2.4p . The 
perturbation equation (14, 2D is globally well-posed in L 1 (M 2 ), and it can be proved that the 
origin w = is a stable equilibrium, but in such a large function space it is impossible to obtain 
a more precise description of the long-time behavior of the solutions. Indeed, it follows from 
the results of [30\ [3Tj that the spectrum of the linearized operator C — a A in L 1 (M 2 ) is exactly 
the full left-half plane {z G C | Re(z) < 0}, for any value of a G M. This means that there exist 
perturbations w which converge to zero at an arbitrarily slow rate. 

For a more precise stability analysis, we have to assume that the perturbations w have a 
faster decay as |£| —> oo than what is strictly needed for integrability. From now on, we write 
w, v instead of w,v, and we consider the perturbation equation (I4.2D in the weighted space 
X = L 2 (M 2 , G~ l d£) equipped with the scalar product 

(w 1 , W2 }= [ G^rwewodc, (4.4) 

and the associated norm ||u>||x = (w, w) 1 ' 2 . Functions in X are locally in L? but have a Gaussian 
decay at infinity; in particular X ^-> L P (R 2 ) for all p G [1,2]. This function space is convenient 
because, as we shall see, the linearized operator C — a A has very nice spectral properties in X. 
However, if we do not not want to restrict ourselves to perturbations with Gaussian decay at 
infinity, it is possible to use the polynomially weighted space l?(ra) = L 2 (M 2 , (1 + |£| 2 ) m d£), 
for some m > 1 [31] . 

The choice of the Gaussian space X is justified by the following result : 

Lemma 4.1 [301 EH [50] 

i) The operator C is selfadjoint in the space X = L 2 (M. 2 , G -1 d£), with compact resolvent and 
purely discrete spectrum : 

°(C) = {-~|n = 0,l,2,...} . (4.5) 

ii) The operator A is skew-symmetric in the same space X : 

(Aioi , w 2 ) + (itfi , Aw 2 ) = , for all wi,w 2 G D(A) C X . (4.6) 

Proof. To prove i), we conjugate C with the square root of the weight G to obtain the operator 



16 ■ 2 ■ (4J » 



L = G- 1 ' 2 C G 1/2 = A- C + I 

which we have to consider as acting on L 2 (M 2 ). Clearly L is the two-dimensional harmonic 
oscillator, which is known to be self-adjoint in L 2 (R 2 ) with compact resolvent and discrete 
spectrum given by (|4.5j) . Returning to the original operator, we obtain the desired conclusions, 
together with the following characterization of the domain of C : 

D(C) = {w G X Aw G X , f • Vu> G A"} . 

The proof of ii) is a direct calculation. By (I4.3P we have A = Ai + A2, where Aiw = v • Vro 
and A 2 w = v-VG = (K*w)- VG. If w 1} w 2 G D(C) C X, then 

(Aiwi,w 2 ) + (wi,A\w 2 } = / G~ 1 (w 2 v ■ Vw\ + w\ v ■ Vw 2 ) d^ 

G~ x v G -V{w lW2 )&i = , 



/ 
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because G l v G is divergence-free. Moreover, since VG = — g£G, we have 

{K 2 wi,w 2 ) + (wi,A 2 w 2 ) = -- / ( (£ • ui)w2 + (£ • v 2 )wi ) d£ 

see also (|2.10|) . Thus (Awi,w 2 ) + (wi,Aw 2 ) = for all wi,w 2 £ D(C) C X. By density, this 
relation holds for all w\,w 2 in the domain of A, which can be characterized precisely [50J. One 
can also show that A is not only skew-symmetric, but also skew-adjoint in X, see |5(Jj . □ 

It is useful to list the first few eigenfunctions of C in X, and to give explicit formulas for the 
corresponding spectral projections. 

1) The first eigenvalue Ao = is simple, and the associated eigenfunction is of course the 
Oseen vortex profile : CG = 0. The corresponding spectral projection Pq : X — > X is given by 
Pqw = G L 2 w d£. We denote by Xq C X the kernel of Pq, namely 



Xn 



\w£X I w &i = 0) = {G} x , (4.8) 

L JR-2 > 



where {G} denotes the set of all w G X which are orthogonal to G in X. Due to the conservation 
of the total circulation, this subspace of X is invariant under the evolution defined by the full 
equation (|4.2p . 

2) The second eigenvalue Ai = —1/2 has multiplicity two, and the associated eigenfunctions 
are the first order derivatives d±G and d 2 G. The spectral projection Pi : X — > X satisfies 
P\w = —diGj R2 ^l^d^ — d 2 G J R2 ?2^d^. Let X\ = ker(Po) H ker(Pi), namely 

X 1 = \weX [ £iiv d£ = fori = l,2\={G;diG;d 2 G} ± . 
L Jr2 y 

Due to the conservation of the first-oder moments of the vorticity, see (|2.8p , the subspace X\ is 
also invariant under the evolution defined by (14. 2p . 

3) The third eigenvalue A2 = — 1 has multiplicity three, and the associated eigenfunctions are 
the second order derivatives AG, {d\ — d 2 )G, and d\d 2 G. Let P 2 be the spectral projection 
associated only with the symmetric eigenfunction AG, i.e. P 2 w = |AG/ r2 |£| 2 wd£. We denote 
X 2 = ker(P ) n ker(Pi) n ker(P 2 ), namely 

X 2 = {wtX 1 I \£\ 2 wdZ = o} = {G;8 1 G;d 2 G;AG} x . 

In view of (12. 9p . the subspace X 2 is also invariant under the evolution defined by (14. 2p . 

More generally, for any n G N, the n eigenvalue A n = —n/2 of C has multiplicity n + 1, 
and the corresponding eigenspace is spanned by the Hermite functions of degree n which can be 
expressed as homogeneous n th order derivatives of the Gaussian profile G. 

We now consider the spectrum of the linearized operator C — aA in the (complexified) 
Hilbert space X, for any fixed a£R. As is easily verified, the operator A is a relatively compact 
perturbation of C in X. Indeed, if we decompose A = Ai + A2 as in the proof of Lemma 14. 11 we 
see that A2 is a compact operator in X, while Ai = v G • V is relatively compact with respect to C 
because this is a first order differential operator whose coefficients decay to zero as |£| — >• 00. By 
a well-known perturbation argument [H] , this implies that C — a A has itself a compact resolvent 
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in X. In particular, its spectrum is a sequence of (complex) eigenvalues {A n (a)} n£ N with finite 
multiplicities, which can accumulate only at infinity. In fact, it is not difficult to verify that 
Re(A„,(a)) — > — oo as n — > oo, for any a G R. 

In general, the eigenvalues of C — ah depend in a nontrivial way on the circulation parameter 
a and are not explicitly known. However, if w G X is a radially symmetric eigenfunction of 
C, then obviously hw = 0, hence w is also an eigenfunction of C — ah for any a G R, and 
the corresponding eigenvalue does not depend on a. In particular, Ao = is an eigenvalue of 
C — ah with eigenfunction G, and the orthogonal subspace Xq = {G} is invariant under the 
linear evolution generated by C — ah. Moreover, if we differentiate the identity v G ■ VG = 
with respect to £i and £2, we see that A(<9jG) = for i = 1, 2. This implies that Ai = —1/2 is a 
double eigenfunction of C — ah for any ct£l, with eigenfunctions d\G and d^G. In addition, 
the orthogonal subspace X\ = {G; d\G; f^G} -1 is invariant under the linear evolution generated 
by C — ah. Finally, since A(AG) = by symmetry, we see that A2 = — 1 is an eigenvalue of 
C — ah for any a£l, and the subspace X2 = {G; d\G; c^G; AG} is also invariant under the 
linearized evolution. 

The following simple but important result shows that the equilibria of fl4.1|) corresponding 
to Oseen vortices are spectrally stable with respect to perturbations in the Gaussian space X. 

Proposition 4.2 [3l] 

For any a£R, the spectrum of the linearized operator C — ah in the space X satisfies 

a{C -ah) C iz g C Re(z) < 0] . (4.9) 

a{C - ah) C iz G C Re(z) < --} in X , 



a{C - ah) C iz G C Re(z) < -l| in X x 



Proof. We already know that the spectrum of C — ah consists entirely of eigenvalues. Assume 
that (C — ah)w = Xw for some A G C and some normalized vector w G D(C) C X. Then, using 
Lemma |4. 11 we find 

Re(A) = Re((£-ah)w,w) = (£w,w) < , (4.11) 

which proves (|4.9|) . Moreover (Cw,w) < —1/2 if w G Xq and {Cw,w) < —1 if w G X\, and 
(I4TT0H follows. □ 

Remark 4.3 // we restrict ourselves to the invariant subspace X2, it is easy to show, as in the 
proof of Proposition ^.^ that all eigenvalues of the linearized operator C — ah satisfy Re(A) < — 1 
i/a/0, but it is difficult to give an optimal upper bound on the real part of the spectrum. 

As an immediate consequence, we can show that the Oseen vortices are linearly stable with 
respect to perturbations in X : 

Corollary 4.4 For alia G R, the semigroup generated by the linearized operator C— ah satisfies 

\\e T{c - aA) \\z^z < e'^ , for all r > , (4.12) 

where /i = if Z = X , fi = 1/2 if Z = Xq, and fj, = 1 if Z = X\ . 
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Proof. We know that the self-adjoint operator C is the generator of an analytic semigroup in 
X, and since A is a relatively compact perturbation of C it is clear that C — a A is also the 
generator of an analytic semigroup, for any a£l. Now, estimate (|4.1ip means that C — ah. is 
m-dissipative in X [41]. By the Lumer-Philips theorem [64], the associated semigroup satisfies 
the bound (|4.12j) with Z = X and /j, = 0. Applying the same argument to the operator C—aA+fj, 
restricted to Xq (if // = 1/2) or X\ (if fi = 1), we obtain the desired result in the other cases. □ 
With some additional work, one can also prove that the Oseen vortices are nonlinearly stable 
with respect to pertubations in X. Here we can restrict ourselves, without loss of generality, to 
perturbations with zero mean; i.e., we can study Eq. (|4.2p in the invariant subspace Xq C X 
defined by (14. 8h . Indeed, adding a perturbation with nonzero mean to the equilibrium aG is 
equivalent to adding a perturbation with zero mean to some modified equilibrium aG, with a 
close to a. The result is : 

Proposition 4.5 [3H [29] 

There exists e > such that, for any a G M., we have the following result. If wq G Xq satisfies 
\\wo\\x < c> then the unique solution of ()4.2|) with initial data wq satisfies 

\\w(t)\\ x < min(l,2e-^ r )||u'o||x , f° r all r>0, (4.13) 

where fi = 1 if wq £ X\ and [i = 1/2 if wq G Xq \ X\ . 

Proof. Arguing as in Section 12. 3| it is not difficult to verify that the perturbation equation 
(|4.2p is locally well-posed in the Gaussian weighted space X. Given wq G Xq with ||u>o||x < e, 
let w G C°([0, T*),X) be the maximal solution of (14. 2D with initial data wq. We know that 
w(t) G Xq for all r G [0,T^), and that w(t) G X\ if wq G X\. To control the time evolution of 
w(t), we use the energy estimate 

5-?-Kr)|& = Wr),£«(r))- / G-^M&rM&r) • Vw(£,r)d£ 

^ dr 7 R2 

= Kr),&(r)) + i / Gr^Xe-vfcrJMCifde, 

4 JR 2 

where in the first equality we used the skew-symmetry of the operator A, and in the second one 
we integrated by parts and used the identity VG _1 = {£,/2)G~ l . 

Fix r > and let /(£) = G _1/2 (0«>(£,i-), for all £ G M 2 . Then \\f\\ L 2 = ||«;(r)||x and, in 
view of P~T[) . 

£(/) := llV/Hi, + ^ ||e/ll| 2 - ^||/||| 2 = -(w(t) , £«;(t)) . 

Since u>(t) G Xo, we have E(f) > /u||/||? 2 where /^ = 1 if wo G Xi and /U = 1/2 if wq G Xq \Xi. 
In particular, there exists C > such that 

£(/) >G(||V/||| 2 + ||C/||| 2 + ||/||| 2 ). 

On the other hand, 

Jr 2 Jr 2 

< \Hr)\\ L ,\\f\\ L ,Uf\\L^ (4-14) 
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By (|2.5p we have \\v\\ L i < C\\w\\ L 4/3 < C||w||x, hence the right-hand side of (|4.14p is bounded 
by C\\w(T)\\xE(f). As a consequence, there exists C\ > such that 

l^\Hr)\\x < -S(/)(l - C!||«;(t)[U) < -»\\w(t)\\ x (1 - Ci||ti;(t)|W , (4.15) 

where the second inequality is valid provided Ci||w(r)||x < 1. 

We now assume that e > is small enough so that C%e < -^. Then the differential inquality 
(J4.15P implies that the norm ||w(r)||x is nonincreasing with time, hence the solution of (|4.2p is 
global (i.e., T* = oo) and satisfies (J4.15P for all r > 0. Integrating (|4.15p . we obtain 

IHr)l|x < . "r°! k ■■ e~>" , forallr>0. 



1-<7i|Kt)|U " 1-Cilkol 



A 



This implies (|4.13p . since Ci||w;o||x < \ by assumption. D 

Proposition 14.51 shows that the steady state w = aG of the rescaled vorticity equation (|4.ip 
is stable with respect to perturbations in the space X, for any value of the circulation parameter 
a £ R. Moreover, we have a lower bound on the size of the local basin of attraction which is 
uniform in a. Here we mean by "local basin of attraction" a neighborhood of the equilibrium aG 
which is small enough so that the time evolution of the perturbations in that neighborhood can 
be controlled by the dissipative properties of the linearized operator, as in the above proof. Of 
course, Theorem l3.1l shows that all trajectories of (j4.1h in L 1 (R 2 ) converge to the line of equilibria 
{aGjagK as r — > oo, hence one can argue that the basin of attraction of these equilibria has, 
actually, infinite size. 

If the initial data are sufficiently localized, Proposition 14.51 allows to improve the conclusion 
of Theorem 13.11 by specifying sharp convergence rates. Assume for instance that the initial 
vorticity uq : R 2 — > R satisfies 

/ / \T\ \ 

I ujq(x) 2 exp(- Ida; < oo , (4-16) 

for some to > 0, and let 7 = L 2 ojq(x) dx. If 7 / 0, we define the center of vorticity as 

xq = — xuJo(x)dx £ R 2 , 

7Jr2 

and we set xq = otherwise. Using the change of variables (|3.14p with xq as above and T = to, 
we transform (|2.2p into the rescaled vorticity equation (|4.ip . and hypothesis (|4.16p means that 
the initial data wq for (|4.1|) belong to the Gaussian space X. Then estimate (|3.18|) (with 
j3 > 1/2) implies that the solution w(t) of (|4.ip is uniformly bounded in X for all r > 0, and 
arguing as in Section [3.31 it is not difficult to show that ||w(t) — aG||x — > as r — > 00, where 
a = 7/V. Thus we can apply Proposition 14.51 when r is sufficiently large, and we conclude that 
||w;(t) — aGr||x = 0(e~^ T ) as r — > 00, for some /x > 1/2. More precisely, when a/0, the choice 
of xq above implies that wq — aG £ X±, hence we can take /J, = 1; otherwise, we take /j, = 1/2. 
Returning to the original variables, we obtain in particular the estimate 

oj(x,t)- , 7 g(^Z^L=) =0(t~ lx ), as t^oo, (4.17) 



v{t + t Q ) \y/v(t + t ) 

which improves (|3.24|) . The decay rate in (|4.17|) is optimal in general, but can sometimes be 
improved by an appropriate choice of to, see [31 (. I29j. We also mention that the convergence 
(|4.17p still holds if, instead of (|4.16p . we assume that the initial vorticity satisfies J R2 ^o(x) 2 (l + 
|x| 2 ) m dx < 00 for some m > 3, see [31]. 
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To conclude our discussion of Proposition 14.51 we would like to emphasize that the stability 
analysis of Oseen vortices lead to conclusions which differ radically from what is known for 
other classical examples in fluid mechanics, such as the Poiseuille flow in a cylindrical pipe or 
the Couette- Taylor flow between two rotating cylinders [23j [78]. Indeed, in the case of Oseen 
vortices, no hydrodynamic instability, neither of spectral nor of pseudospectral nature, develops 
as the circulation parameter a = j/v (which plays the role of the Reynolds number) is increased. 
On the contrary, it is possible to show that a fast rotation has even a stabilizing effect on the 
vortex, as we shall see below. 

4.2 Improved Stability Estimates for Rapidly Rotating Vortices 

For a more detailed study of the linearized operator C — oA, it is useful to observe that both 
operators C and A are invariant under rotations about the origin in R 2 . This symmetry is fully 
exploited if we introduce polar coordinates (r, 8) in R 2 and expand the vorticity distribution in 
Fourier series with respect to the angular variable 6. Let us decompose 

X = X n , (4.18) 

nGN 

where X n denotes the subspace of all w £ X such that 

w(r cos 8, r sin 8) = a(r) cos(n0) + b(r) sin(n#) , r > , 8 £ [0, 2tt] , 

for some radial functions a, b : R+ — > R. The invariance under rotations implies the following 
result, which can be obtained by a direct calculation : 



Lemma 4.6 [31] 

For each n E N the subspace X n is invariant under the action of both linear operators C and A. 
The restriction C n of £ to X n is the one- dimensional operator 

The restriction A n of A to X n satisfies A n uj = in(<poj — g£l n ), where 

' :i-e" r2 / 4 ), g(r) = ^e-r 2 /\ r>0, (4.19) 



2irr 2 ' Att 

and (for n > 1) Q n is the regular solution of —Q'n Q' n + \Q n = u, namely 

Q n (r) = -!- ( f (-Yr'uj{r')dr' + f (^Yr'uo(r')dA , r>0. (4.20) 

In particular, Lemma 14.61 implies that the subspace Xq C X, which consists of all radially 
symmetric functions, is entirely contained in the kernel of the operator A. As was already 
observed, the first order derivatives diG,d2G also belong to ker(A), and the following result 
shows that these are the only non-radially symmetric functions in the kernel of A : 

Lemma 4.7 [50] ker(A) = X {MxG + /3 2 d 2 G | Pi, fa € R}. 

Proof. We already know that Xq C ker(A). Suppose now that w G ker(A) DX n for some n > 1. 
Without loss of generality, we can assume that w{r cos8,r sin#) = uj{r) cos(n8) or ui(r)sin(n8) 
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for some uj : R+ —> R, and using the notations (|4.19p . (|4.20p we have (jxv — g£l n = 0. This means 
that $7„ satisfies the differential equation 

" Kir) - Itf B (r) + (^ - f|})^(r) = , r > . (4.21) 

Now it is easy to verify that r 2 g(r)/cj)(r) < 4 for all r > 0. Thus Eq. (|4.2ip satisfies the maximum 
principle if n > 2, hence has no nontrivial solution vanishing at the origin and at infinity. This 
shows that ker(A) n X n = {0} for n > 2. When n = 1, Eq. (|4.2ip has a (unique) nontrivial 
solution given by fii(r) = r(f>(r), which satisfies (|4.20p with uj{r) = rg(r). This shows that 
ker(A) n X\ is spanned by d\G and c^G. D 

In the rest of this section, to investigate the effect of the rotation on the stability of Oseen's 
vortex, we restrict ourselves to the orthogonal complement of ker(A) in X. This subspace, 
denoted by ker(A)^, is invariant under the evolution defined by the linearized equation d T w = 
(C — aA)w, but not under the full equation (I4.2p . Let C± and Aj_ denote the restrictions of 
the operators C and A to the invariant subspace ker(A)- 1 -. For any aelwe define the spectral 
lower bound 

S(a) = infJRe(z) z G a(-C± + aA±)\ , 

and the pseudospectral bound 

tf(a) = (sup||OCj_-aA._L-tA)- 1 ||x_ K x 



Note that, in the definition of S(a), we have changed the sign of the linearized operator, in 
order to get a positive quantity. 

Lemma 4.8 For any a£l one has £(a) > \P(a) > 1. 

Proof. Fix a G M. Since ker(A)- 1 C X h it follows from (j4~10j) that £(a) > 1. Moreover, if 
(C — aA + z)w = for some z £ C and some normalized w G ker(A)^, then (>C— QA + iIm(z))7i; = 
— Re(z)w, so that 

Re(z) > \\(C ± -aA ± + ilm(z))- 1 \\- 1 > *(a) . 

This proves that S(a) > ^/(a). On the other, the proof of Proposition 14.21 shows that the 
operator C± — aA± + 1 is m-dissipative. This in particular implies that || (C± — aA± — iX)^ 1 \\ < 1 
for all Ael, hence *(a) > 1. □ 

Numerical observations due to Prochazka and Pullin |65| indicate that X(a) = Oflal 1 ' 2 ) 
as \a\ — > oo. This has not been proved yet, but a recent result by Maekawa [50] shows that 
£(a) —> +oo as \a\ —> oo. Although the proof is not constructive and does not give any 
information on the asymptotic behavior of X(a), this is the first rigorous result which shows 
that a fast rotation stabilizes Oseen's vortex, in the sense that the decay rate of the perturbations 
in ker(A)^ becomes arbitrarily large as \a\ — > oo. In fact, the argument of [50] can be easily 
modified to yield the stronger result \I/(a) — > oo as \a\ — > oo. This indicates that the size of the 
local basin of attraction of Oseen's vortex aG becomes arbitrarily large (in the Gaussian space 
X) when \a\ —> oo. 

Explicit bounds on £(a) or ^(a) are more difficult to obtain, and were first established 
in [26] for a model problem in one space dimension, and in [21] for a simplified version of the 
linearized operator C — aA where the nonlocal term A2 in A is omitted. In both cases one gets 
ty(a) = 0(\a\ ' 3 ) as \a\ — > 00, and for the one-dimensional model we also have X(a) = C(|a| ' ). 
These partial results suggest that S(a) 3> ^(a) when \a\ is large, a property which reflects the 

30 



fact that the linearized operator is highly non-selfadjoint in the large circulation regime. Indeed, 
it is easily verified that the spectral and pseudospectral bounds always coincide for selfadjoint 
operators. 

Although the analysis of the linearized operator C — aA is still the object of ongoing research, 
we formulate the following conjecture which should be proved soon. 

Conjecture 4.9 There exist constants 5 > 7 > such that ty(a) = C(|a| 7 ) and S(a) = C(|a| ) 
as \a\ — > 00. 

On the basis of the analysis of simpler models [26], we further conjecture that 7 = 1/3 and 
5 = 1/2, but these optimal values may be difficult to reach. To conclude this discussion, we 
mention that, if the linearized operator C — a A is restricted to the invariant subspace 

© X n C ker(A) ± , 

n>N 

for some sufficiently large integer N > 2, then it is possible to show that the pseudospectral 
bound in that subspace satisfies \&Ar(a) = Cda) 1 ' 3 ) as \a\ —> 00, see J22]. 

5 Interaction of Vortices in Weakly Viscous Flows 

Theorem 12.81 shows that the two-dimensional vorticity equation (12. 2\i is globally well-posed in 
the space of finite measures, for any value of the viscosity parameter v. It is therefore natural 
to investigate how the solutions of (12.2p behave in the vanishing viscosity limit. This question is 
very difficult, because the formal limit of fl2.2f) as v — > is the two-dimensional inviscid vorticity 
equation 

d t co(x, t) + u(x, t) ■ Vuj(x, t) = , (5.1) 

which is certainly not well-posed in the space A^(M 2 ). In fact, a celebrated result of Yudovich [81] 
shows that (15.1J) has a unique global solution if the initial vorticity \i belongs to L 1 (R 2 )nL oc (M 2 ), 
see also [821 [80]. If ^ G M(M. 2 ) n i? _1 (lR 2 ) and if the singular part of fi has a definite sign, then 
(|5.ip has at least one global weak solution |20|. I5T] , but this result does not cover the case where 
the initial data contain point vortices. 

A careful treatment of the inviscid limit is beyond the scope of these notes, but the following 
brief discussion will be useful to put our results into perspective. As a general rule, it is known 
that the Navier-Stokes equation is nicely approximated by the Euler equation in the vanishing 
viscosity limit if we restrict ourselves to smooth solutions in a domain without boundary [74 1 142 1 
[2]. Thus the main difficulties come either from the presence of boundaries or from singularities 
in the initial data. Since we chose to work in the whole plane M 2 , we do not mention here the 
(hard) problems related to boundary layers, and refer to [69j [70J [35] for a discussion of that 
point. But even in the absence of boundaries, the inviscid limit may be nontrivial to understand 
if one considers non-smooth initial data. In the case of vortex patches, the situation is well 
understood : the Euler flow is well defined due to Yudovich's theorem, and approximates the 
Navier-Stokes flow in every reasonable sense asz^^O, see [HI [151 HH HH HH [H EH EH1 EH1 173] . 
The situation is much less clear for vortex sheets : although an inviscid solution can still be 
constructed, the Euler flow is unstable due to the discontinuity of the velocity field (Kelvin- 
Helmholtz instability), and the vanishing viscosity limit is as difficult to study as in the case of 
Prandtl boundary layers [TJ. Finally, the case of point vortices is somewhat paradoxical : one the 
one hand, these are the most singular initial data for which we can solve the 2D Navier-Stokes 
equation, and strictly speaking the Euler flow is not even defined for such initial data. On the 
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other hand, we have a nice substitute for the Euler flow in that particular case, namely the 
Helmholtz-Kirchhoff point vortex system, which does not exhibit any dynamical instability. As 
a result, the inviscid limit can be rigorously studied in the case of point vortices. Early results in 
this direction were obtained by Marchioro [MIES], and will be described more precisely below. 

5.1 The Viscous N- Vortex Solution 

From now on, we fix an integer N > 1, and we assume that the initial vorticity is a collection of 
iV point vortices characterized by their (pairwise distinct) positions x±, . . . , xn £ K 2 and their 
(nonzero) circulations 71 , . . . , 7^ 6t: 

N 

H = 5^7t<*s* • (5.2) 

i=i 

Given a viscosity v > 0, we denote by uj y (x,t), u u (x,t) the unique solution of the vorticity 
equation (j2.2[) with initial data /U. Existence of such a solution was first proved in [4], and 
uniqueness is guaranteed by Theorem 12.81 Note that Theorem 12.51 applies if ^ | -y^ | < Cqv, but 
such an assumption is totally inappropriate here because we want to study the limit v — > for 
fixed initial data. 

If N = 1, the solution of (12. 2|) is just the Oseen vortex with circulation 71 centered at the 
point x\ : 

v, ,n 7i r,( x-x\ \ v, ,n 7i G (x-x\ 

to (x,t) = — Gl — -==-) , u (x,t) = —j=v ' 



vt V y/vt / yvt V yvt 

When N > 2, the viscous N- vortex solution is no longer explicit, but as long as the diffusive 
length \ful is small compared to the distance between the vortex centers we expect that uj u (x,t) 
can be approximated by a superposition of Oseen vortices : 

«"M>*x:^(^f), „ M) ,£^ ( ^f), ,,3, 

i=i v i=i v v 

where z±(t), . . . , z^{t) denote the positions of the vortex centers at time t. In the same regime, 
we also expect that the positions Zi(t) will be determined by the Helmholtz-Kirchhoff system 

[SEE]: 

<w = ■kX-« M)~-tw ' m = x - (5 ' 4> 

Indeed, we know from the work of Marchioro and Pulvirenti that the system of ordinary dif- 
ferential equations (|5.4p approximates rigorously the motion of localized vortex patches in two- 
dimensional inviscid [53], [56j [57] or slightly viscous [5J1 [55] fluids. 

The point vortex system (|5.4p is globally well-posed for all initial data if iV = 2 or if all 
circulations 7^ have the same sign. In the general case, however, vortex collisions can occur 
in finite time for some exceptional initial configurations |57l 62J. To eliminate this potential 
problem, we assume in what follows that the solution {zi(t), . . . ,ZN(t)} of (|5.4h is defined on 
some time interval [0, T], and we denote by d the minimal distance between the vortex centers : 

d = min min Izdt) — Zj(t)\ > . (5-5) 

te[o,T] i+j 

We also introduce the turnover time To = d 2 /\j\, where I7I = I71I + • • • + |7at|, which is a natural 
time scale for the inviscid dynamics described by (|5.4p . For instance, for a pair of vortices with 
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the same circulation 7 separated by a distance d, one can check that the rotation period of each 
vortex around the midpoint is 4tt 2 Tq. 

With these notations, we can formulate now the main result of this section : 



Theorem 5.1 [28] 

Assume that the point vortex system (|5.4p is well-posed on the time interval [0, T]. Then the 
unique solution of the two-dimensional vorticity equation (|2.2p with initial data fi = Yli=i 7* $b* 
satisfies 

where {zi(t), . . . ,zj^(t)} is the solution of (|5.4|) and K is a dimensionless constant depending 
only on the ratio T/Tq. 

Theorem 15.11 describes the viscous iV- vortex solution to leading order in our expansion pa- 
rameter \puijd, which is the ratio of the typical size of the vortex cores to the minimal distance 
between the vortex centers. In other words, the approximation (|5.3|) is accurate whenever the 
vortices are widely separated compared to their size. Another dimensionless quantity that is 
present in our analysis is the ratio of the observation time T to the turnover time To . It is im- 
portant to remark that no smallness assumption is made on this ratio in Theorem 15.11 although 
the constant K in (j5.6|) becomes very large if T 3> To . 

Since each Oseen vortex converges weakly to a Dirac mass in the inviscid limit, an immediate 
consequence of Theorem 15.11 is : 

Corollary 5.2 Under the assumptions of Theorem \5. 11 the viscous N -vortex solution to v {x,t) 
satisfies 



N 

1/ / ■ \ 

1=1 
where {zi(t), . . . , zjy(t)} is the solution of (|5,4p . 



(.,*) -* t>i<S*(t), forallte[0,T], (5.7) 



While weaker than Theorem 15. 1( this corollary provides a rigorous derivation of the point 
vortex system (J5.4J) using the Navier-Stokes equation in the inviscid limit. In contrast, the classi- 
cal approach of Marchioro and Pulvirenti [56^157] allows to justify (I5.4p within the framework of 
Euler's equation, but requires an approximation argument which is not needed in Corollary 15.21 

5.2 Decomposition into Vorticity Profiles 

In the rest of this section, we give a sketch of the proof of Theorem 15.11 Our starting point is 
a natural decomposition of the viscous iV-vortex solution, which follows immediately from the 
representation formula (|3.10p . Since the initial measure \i is a sum of Dirac masses, we have by 
(|3lT)jl 

N N 

u u (x,t) = £>?(*, t) , u"(x,t) = £>r(M) , (5-8) 

i=l i=l 

where 

u\{x,t) = ^K(x,t; Xi ,0) , xel 2 , 0<t<T, i = l,...,N, 
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and u\(x,t) is the velocity field corresponding to u) v {x,t) via the Biot-Savart law (|2.4p . Here 
r^(x, t; y, s) is the fundamental solution of the convection-diffusion equation (13. 4p with advection 
field U(x,t) = u u (x,t). In view of (J3.6J) . we have the bound 

Wix,t)\ <K^ exp(-^^|£) , (5.9) 

for all i G {1, . . . ,N}, all a; G M 2 , and all t G (0,T], where f3 G (0, 1) and the dimensionless 
constant K\ > depends only on the total circulation I7I and on the viscosity v. Estimate (|5.9p 
gives a precise information on the iV-vortex solution when v > is fixed and t — > 0, but cannot 
be used to control u^(x, t) when t > is fixed and v — > 0, because the constant K\ blows up in 
that regime. 

To analyze more precisely each term in the decomposition (|5.8p we introduce, for each i G 
{1, . . . , N}, the self-similar variable 

X-Zi(t) 

? = - 



vt 

where {zi(t},..., Zjv(i)} is the solution of (|5.4p . As in (]3.14p . we define rescaled vorticities 
w%(£,t) G R and velocities <(£,£) G M 2 by setting 



Our goal is to describe the vorticity profiles w" (£, i) as precisely as possible, in the regime where 
\fvi < d. According to (pTfijh we have <(^i) = G(^) + 0(vt/d 2 ) for each i G {l,...,iV}, 
where G is the profile of Oseen's vortex. However, it is very important to realize that this 
leading order approximation is not sufficient to control the inviscid limit of the viscous iV- 
vortex solution. The reason is that the profile G is radially symmetric, and thus does not take 
into account the deformations of the vortices due to mutual interactions. 

If we go back to the explicit formula (|3.ip . we see that the velocity field of Oseen's vortex is 
very large near the center if the viscosity u is small, with a maximal angular speed of the order 
of \"f\/{vt). As long as the vortex stays isolated, this large velocity does not affect the vorticity 
distribution, which is perfectly radially symmetric. Now, if the same vortex is advected by a 
non-homogenous external field, which in our case will be the velocity field produced by the other 
N— 1 vortices, its vorticity profile gets deformed under the external strain and the vortex starts 
feeling the effect of its own velocity field. As is easily verified, this self-interaction effect is very 
strong if the circulation Reynolds number I7I/V is large, even for a moderate deformation of the 
vortex core. This gives a clear indication that self-interaction effects must be taken into account 
if one wants to control the vanishing viscosity limit in presence of point vortices. 

The deformation of a rapidly rotating Oseen vortex advected by an external velocity field 
can be computed using the properties of the linearized operator C — aA, which was defined 
in Section [U Remarkably enough, one finds that Oseen's vortex adapts its shape in such a 
way that the self- interaction counterbalances the effect of the external field |77j 176] . except 
for radially symmetric corrections and for a rigid translation. This is related to Lemma 14.71 
which shows that the kernel of the operator A consists of all radially symmetric functions and 
of the two-dimensional subspace spanned by d\G, c^G. This fundamental fact explains why 
one can observe, in turbulent two-dimensional flows, stable asymmetric vortices which in a first 
approximation are simply advected by the main stream [5]. 

In Section 15.31 we implement this idea by constructing a higher order approximation of 
the viscous iV-vortex solution which takes into account the self-interaction effects. For each 
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i £ {1, . . . , iV} and all t £ [0, T], the approximated vortex profile will be of the form 

< PP (£,i) = G(0 + Q){m,t) + F?(ti,t)} , £eM 2 , (5.11) 

where Fi(£,t) is a radially symmetric function of £ which we will not describe precisely because 
it only represents a small perturbation of the Gaussian profile G(£). The important correction 
at this order is the nonsymmetric function i^"(£, t) which is given by 

where we denote 2«(£) = Zi(t) — Zj(t). Here w : (0,oo) — > R is a smooth, positive function 
satisfying w(r) ~ Cir 2 as r — >• and ui(r) ~ C2r 4 e _r ' 4 as r — >• oo for some C\, Ci > 0. 

The right-hand side of (|5.1ip is the beginning of an asymptotic expansion of the rescaled 
vortex patch w"{^,t) in powers of the non-dimensional parameter x/ui/d, which is the ratio of 
the size of the vortex cores to the minimal distance between the centers. Each term in this 
expansion can in turn be developed in powers of the inverse circulation Reynolds number y/\"f\. 
The most important effect is due to the nonsymmetric term F"(£,t), which describes to leading 
order the deformation of the i th vortex due to the influence of the other vortices. Keeping only 
that term and using polar coordinates £ = (r cos#,rsin#), we can rewrite (15. lip in the following 
simplified form 



(«,«) = sM + ^Ef i^cos^-M*))) + ... , P.13) 



li \ z ij{t)\ 



where g(\£\) = G(£) and Oij(t) is the argument of the planar vector %■(£) = Zi(t) — Zj(t). This 
formula allows us to compute the principal axes and the eccentricities of the vorticity contours, 
which are elliptical at this level of approximation. 

To formulate our result in its final form, we introduce a function space that is appropriate 
to control the vorticity profiles. Given a small /3 £ (0, 1), which will be specified later, we define 
the weighted L 2 space Xp equipped with the norm 

\H\x g = (f \w(Z)\ 2 e^\/UA l/2 . (5.14) 



In particular, we have Xp <-t L 1 (M 2 ) for any /3 > 0. As is explained in [28], it would be more 
natural to use here the Gaussian space X defined in (14. 4\\ . but for technical reasons it appears 
necessary to replace the Gaussian weight by an exponential one, at least if the ratio T/Tq is 
large. 



Proposition 5.3 Assume that the point vortex system (|5.4p is well-posed on the time interval 
[0, T\, and letu u (x,t) be the solution of (|2.2p with initial data (15. 2p . Then there exist positive 
constants (3 and Ki, depending only on the ratio T/Tq, such that, ifuj u (x,t) is decomposed as 
in (|5.8j) . then the rescaled profiles w^(^,t) defined by (I5.10p satisfy 

.maxJKM) -uF*{;t)\\x, < K iQ) 3/2 . (5-15) 

for all t £ (0,T], provided uT/d 2 < K^ 1 ■ 
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As is clear from (|5.15p . (|5.1ip . and (|5,10p . Proposition 15.31 implies immediately Theorem 15.11 
Note that the error term on the right-hand side of (J5.15J) is smaller than in (15. 6p , and in particular 
smaller than the first order corrections to the Gaussian profile in ()5.1ip . This means that 
the deformations of the interacting vortices are indeed described, to leading order, by (|5.13p . 
According to that formula, each vortex adapts its shape instantaneously to the relative positions 
of the other vortices, without oscillations or inertia. By this we mean that, for each t G (0, T], 
the angular factor cos(2(# — Oij)) in ()5.13p . which gives the leading order deformation up to a 
time-dependent prefactor vt/\zij\ 2 , is entirely determined by the instantaneous positions of the 
vortex centers. In contrast, it is shown in [28j that the first order radially symmetric corrections 
Fi(£,t) do not only depend on the instantaneous vortex positions, but on the whole history of 
the system. 

It is instructive at this point to recall what is known in the more general situation where 
the initial data are not point vortices, but finite size vortex patches. For simplicity, we consider 
the case of two identical, radially symmetric vortex patches of size R > initially located at 
points xi,X2 £ M 2 with \x± — X2I 3> R. Under the time evolution defined by (J2.2J) . the centers 
of the patches will rotate with approximatively constant angular velocity around the mid-point 
(xi + X2V2. To obtain a more precise description of the solution, we go to a rotating frame 
where the centers of the patches stay fixed, and we concentrate on the deformations of the vortex 
cores. What is observed in numerical simulations [45j is that the interaction begins with a fast 
relaxation process, during which each vortex adapts its shape to the velocity field generated by 
the other vortex. This first step depends on the details of the initial data, and is characterized 
by temporal oscillations of the vortex cores which disappear on a non-viscous time scale. In 
a second step, the vortices relax to a Gaussian-like profile at a diffusive rate, and the system 
reaches a "metastable state" which is independent of the initial data, and will persist until two 
vortices get sufficiently close to start a merging process. In this metastable regime, the vortex 
centers move in the plane according to the Helmholtz-Kirchhoff dynamics, and the vortex profiles 
are uniquely determined, up to a scaling factor, by the relative positions of the centers. This is 
exactly the situation described by Proposition 15.31 In fact, if we start with point vortices, the 
system is immediately in the metastable state which, for more general initial data, is reached 
only after the transients steps described above. In this sense, point vortices can be considered 
as well-prepared initial data for the vortex interaction problem. 

5.3 Perturbation Expansion and Error Estimates 

In this final section, we briefly indicate how to construct an asymptotic expansion of the viscous 
N- vortex solution and to control the error terms. We first write the evolution equation satisfied 
by the rescaled vorticity profiles defined in ()5.10p . Since dtoj" +u v - Vaj" = v/S.ui" for i = 1, . . . , TV, 
where u v = X^'=i n ^> we obtain 

N 

" ■ ■ ( / i \ . / 

«) 'V«f(£,t) = (£<)«,() 



««+ lM« + f ■')"# 



where C is the linear operator (|4.ip . The left-hand side is clearly singular in the limit v — > 0, 
but this difficulty can be partially eliminated by an appropriate choice of the speeds z[(t). To 
this purpose, we set 



N a\ 



7=1 



(5.16) 
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where Zij(t) = Zi(t) — Zj(t) and v G is the velocity profile defined in (|3.2p . This is a viscous ap- 
proximation of the original point vortex system (|5.4p . which has its own interest and was recently 
studied in [6TJUD]. It is shown in |28j Lemma 2] that the solutions of (j5.16|) stay extremely close 
to those of (15. 4p in our perturbative regime where the vortices are widely separated from each 
other. In what follows, we thus make no distinction between the solutions of (|5.4p and (|5.16p . 

As was already mentioned, we expect the Gaussian function G to represent the leading order 
approximation of each vorticity profile. Substituting G for w"(£,t) in the evolution equation 
above, and using (I5.16J) together with dtG = CG = 0, we obtain the following expression for the 
residuum of this naive approximation : 

^,^E^ + ^)-^)}.v <o. 

This expression is not singular in the limit v — >• 0, and a direct calculation (see [28^ Proposition 
1]) yields the asymptotic expansion 

Rfi^t) = ^{A t (C,t)+Q) 1/2 B l ^t) + Q)c^,t) + Rf > \^ , (5.17) 

for all £ G R 2 and all t 6 (0,T], where 

MU) ~ ^\- % MOP G ® ' 



<*«,.) = - E S K ' I,J ? )> l;' i(t>1) fat ■ *M) 3 - If I' N Ml') o«) . 



; .. 7i M*)l e 
Moreover, given any (5 < 1, the last term i-Q in (|5.17p can be estimated as follows : 

\Rf\tM < cg) 3/2 e-^i 2 / 4 , £ e M 2 , o < t < T . 

The formula (J5.17P shows in particular that i?| (£,£) = 0(1) as f — > 0. Thus, if we decom- 
pose u^(£,£) = G(£) + Wi(t;,t), the equation for u)j(£,£) will contain a source term of size 0(1) 
as i/ — )■ 0, and we therefore expect that the remainder Wi(£,t) itself will be of size 0(1) after a 
short time. But, as is easily verified, the equation for wij(£, t) contains nonlinear terms involving 
negative powers of v, and such terms cannot be controlled in the vanishing viscosity limit if 
Wi(£,t) is 0(1). This is the reason why it is necessary to construct a more precise approximate 
solution of our vorticity profiles, in order to "desingularize" the equation for the remainder. 

To this end, we first consider approximate solutions of the form 

< pp (£,*) = g(o + (^)*i(e,o , < pp (e,t) = v G (o + i^y^t) > 

where i*i(£, t), . . . , F n (£, t) are smooth vorticity profiles to be determined, and the corresponding 
velocity profiles v Fl (£,, t), . . . , v Fn (£,, t) are obtained via the Biot-Savart law (|2.4p . The residuum 
of this improved approximation can be computed [28J , and has the form 



*!%*) = Itetj+A^til+o'^ 
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where A is the linear operator defined in (|4.3p . To minimize the residuum, we want to choose Fi 
so as to satisfy the linear equation AFi + Ai = 0. This is indeed possible, because A{ belongs to 
the Schwartz class 5(M 2 ) and to the subspace Xi C X, see (|4.18|) . and it is shown in [60l 1321 128] 
that 5(M 2 ) n X2 C ran(A). The equation AFi + A4 = has therefore a unique solution in X2, 
which coincides with the leading term in the expression f|5.12j) of i 7 ^ (£,£). Note however that 
the equation above determines F{ only up to an element of ker(A). As a matter of fact, in (|5.11|) . 
a radially symmetric correction Fj(£, £) is added to F"{^,t) in order to compensate for higher 
order terms (which will not be considered here). In any case, we can choose Fi in such a way 
that 2?j (1) (C,t) = 0{^7t/d). 

Following the same procedure, we next construct a higher order approximate solution of the 
form 



app 



(>«,«) +(|) 3/ W)+ (IP 



< pp (e,t) = g(o + (^)^(e,t) + (^) 3 W,*) + (^) 2 ^(e,*) , (5.18) 



(£,*) = ^0 + (£)**&*) + (£) '>& *) + (£) W,t) 



Again, the vorticity profiles F$, Hi, and Ki are obtained by solving under-determined linear 
equations involving the operator A. More precisely, we have AFi + Ai = and AHi + E>i = 0, 
but the equation for Ki is more complicated. The residuum of that approximation now satisfies 

\R?\^)\ <cg) 3/2 e" W4 , ^f 2 , 0<i<T. (5.19) 

Finally, once the approximate solution (|5.18p has been constructed, we decompose the vor- 
ticity profiles as follows : 

<(£,i) = wf v {i,t) + g) w^t) , <(£,*) = < pp (£,i) + g) *(£,*) , (5.20) 

and we consider the evolution system satisfied by the remainder w)j(£, t), £j (£,£). This system 
reads 

+ ^(% app (e,i) • v^(e,t) + «i(e,t) • v< pp (e,t)) 



j¥* 



AT 



Jjt „ / ^y(t) 



+ E|^(^ + T^' t )- V ^ app ^ f) (5 - 21) 



+ Ef : ^^ + v^' t )- v ^ (e ' t)+ ^ (e ' t) = ' 



3=1 

and using (|5.19p . (|5.20p we see that the residuum Ri is 0(y/ui/d) as v — >• 0. Also, due to the 
ansatz (|5.20p . the nonlinear terms in (|5.2ip are now regular in the limit v — )■ 0. The linear terms 
still contain negative powers of u, but this difficulty can be avoided by using appropriate norms. 
In a much simpler setting, the same idea was already present in the proof of Proposition 14.51 
where the skew-symmetry of the operator A in the Gaussian space X was used to obtain a 
stability estimate that was uniform in the circulation parameter a. 
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To control the solution of (|5.2ip . we introduce a weighted energy functional of the form 

N 

If the observation time T > is small with respect to the turnover time To = d 2 /\^/\, we can 
take Pi(£, t) = p a (t)(0 for i = 1, . . . , N, where a(t) = d/(3yut) and 



p«(0 



el«l / 4 if 1^1 < a , 
e a2 / 4 if a < |f | < ifa 
e |^/(4^ 2 ) if \S\>Ka, 



for some K > 1. We then have e ^ 2 ^ AK ^ < Pi{£,t) < e l^l 2 / 4 for all £ G M 2 and t G (0,T). With 
this choice, we obtain from (|5,2ip a differential inequality for the weighted energy E(t), which 
can be integrated using Gronwall's lemma and yields the bound : 



/ 



<^{\w 1 (^t)\ 2 + --- + \wN^,t)\ 2 )d^ < E(t) <C^ 



This concludes the proof of Proposition 15.31 if T <C To. In the general case, one has to introduce 
more complicated weights, which can be constructed using the same procedure as the approx- 
imate solution itself. These weights satisfy e^'' 4 < Pj(£, t) < e'^' ' 4 , for some small /3 > 
depending only on T/Tq . We thus obtain the weaker estimate : 



/ e ^ i (\w 1 &t)\ 2 + --- + \w N ({;,t)\ 2 )d£ < E(t) < C 



vt 



which still implies the desired conclusion. 
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